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DESIGN CONSIDERATIONS IN THE FIRST STAGE 
OF TRANSISTOR RECEIVERS 

By 

Larry A. Freedman 

RCA Laboratories, 
Princeton, N". J. 

Summary — This paper presents a disc2issw7i of noise performance of 
transistor r-f stages utilizing capacitive antennas and of transistor mixer 
stages utilizing loop antenna^. Examples of the noise performance to be 
expected with each type of antenna are inclnded. Comparisons are drawn 
between transistor stages and corresponding tube stages. Consideration is 
given to design compromises between image rejection and insertion loss for 
an r-f stage employing both tuned input and interstage transformers. The 
procedure for transformer design for optimum iiisertion loss-image rejection 
performance is outlined. 



Noise Performance 

Capacitive Antenna 

ATYPICAL transistor r-f stage is shown in Figure la. A capaci- 
tive-type antenna is assumed (as for example the rod antenna 
-of an automobile receiver). Figure lb shows an equivalent 
circuit from the antenna to the r-f stage input where the antenna has 
been replaced by a voltage source, y i, and an internal impedance con- 
sisting of a capacitance C^. The transistor is represented by its input 
resistance, R^, and the transformer tuned impedance by R^. The signal 
power delivered to the r-f stage input by an unmodulated signal is 
given by 



r 



R 

P 



n^-R, + R, _ 



1 

— , (1) 

R; 



where n is the transformer turns ratio, NJN^, as indicated in Figure 1 
(essentially unity coupling is assumed), Ct is the sum of the antenna 
capacitance and the total shunt capacitance across the antenna trans- 
former, and Qoi is the unloaded Q of the antenna transformer. This 
expression is derived in Appendix I. 

The equivalent thermal noise generator referred to the secondary 
of the equivalent circuit of Figure lb is shown in Figure Ic. The 
antenna radiation resistance is negligible compared to R^ and is ignored. 
The thermal noise power delivered to the transistor is 
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(2) 



where k is Boltzmann*s constant, T is absolute temperature, and A/ is 
the noise bandwidth. Since the r-f transistor noise factor, F, may be 

AMTENNA TRANSFORMER 
Np » PRIMARY TURNS 
Ng » SECONDARY TURNS 

INTERSTAGE TRANSFORMER 
My- TOTAL PRIMARY TURNS 
Hl Mc- primary turns AT COaECTOR 
Mj- SECONDARY TURNS 

MIXER INPUT 
RESISTANCE 

o) TRANSISTOR RF STAGE 




Np: Ns 



b» EQUIVALENT INPUT CIRCUIT 



e *N/4KT(n^R„) Af 



c) EQUIVALENT THERMAL NOISE INPUT CIRCUIT 



Hh- 



d) TUBE INPUT CIRCUIT 

Fig. 1 — R-F stage utilizing capacitive antenna. 

expressed as Py P-rih where Py is the total noise power referred to the 
transistor input, 

l\^FPru (3) 



(assuming that noise from the second and succeeding stages is small 
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and rnay be neg'lected). The signal-to-noise ratio is obtained from 
Equations (1), (2), and (3) and is given by 



I\r AkT Af F 



(4) 



For a receiver employing a linear detector the signal-to-noise ratio at 
the output of the audio amplifiei', (J^R/pN)Ay is related to the r-f 
signal-to-noise ratio by 

P. \ Ps A/ 

, (5) 

I, F^y Af^ 

where 7)1 is the per cent modulation and A/4 is the over-all noise band- 
width of the receiver. (For the usual case where the receiver audio 
bandwidth is less than half of the bandwidth to the detector input, the 
over-all noise bandwidth is twice the audio noise bandwidth.) For a 
30 per cent modulated signal and an output signal-to-noise ratio of 20 
decibels, (P.s^-PaO.^ ~ 100; the antenna voltage required for this signal- 
to-noise ratio is obtained by combining Equations (4) and (5) and is 
given by 



lOV'^kTR^ Af^F 

V,^ (6) 

Ca 

0.3- Qoi 

This expression is retained in this form to facilitate comparison with 
the tube case discussed below (the mutual dependence of R^, Cj., and 
Qoi should be noted, i.e., Qqi =3i oi C^j. R^) . Note that Equation (6) does 
not appear to be directly dependent on the turns ratio or operating Q 
of the antenna transformer. The transistor noise factor, however, is 
dependent on the turns ratio to the extent that the turns ratio deter- 
mines the driving source impedance. This noise factor dependence is 
quite small in the vicinity of the matched source impedance. For 
example, a four-to-one mismatch in either direction causes an increase 
in F of about 2 decibels* corresponding to an increase of about 25 per 
cent in the signal required for 20 decibels signal-to-noise ratio. 

A comparison of the noise performance of a tube and transistor r-f 



* Determined experimentally. 
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stage may be made by comparing Equation (6) with a like expression 
for a tube stage. Figure Id shows a typical antenna connection for a 
tube r-f stage where the primary of the antenna transformer of Figure 
lb is connected directly to the grid. The tube noise factor, F^, may 
be expressed as 

= (7) 

where R^.^^ is the equivalent noise resistance of the tube. For typical 
remote-cutoff pentodes, the equivalent noise resistance is of the order 
of 3,000 to 5,000 ohms. For this type of connection at broadcast fre- 
quencies is usually large compared to 5,000 ohms and Fj^ is very 
nearly unity. The signal required for 20 decibels signal-to-noise ratio 
is then 



~. (8) 

Ca 

0.3 Ooi 

A comparison of Equations (6) and (8) shows that the tube r-f stage 
will be times better than a transistor stage. 

Example of Noise Performance Calculation 

As an example, consider the r-f stage of a receiver for which the 
input circuit parameters are C,i ~ 30 micromicrofarads, Ct = 100 
micromicrofarads, and Qqi = 70. The noise factor of the r-f transistor 
is 6 decibels and the over-all noise bandwidth is 3 kilocycles. Substi- 
tuting the pertinent values in Equation (6), 



lOVl.G X 10--^' X 112 X 10'^ X 3 X 10*^ X 4 

T^i = = 7.37 microvolts 

30 

0.3 X X 70 

100 

for a 20-decibel signal-to-noise ratio at 1 megacycle. 
Loop Antenna 

A mixer or converter first stage is frequently used in receivers 
employing a loop antenna. A transistor mixer first stage is shown in 
Figure 2a. The antenna circuit typically consists of a coil wound on 
a ferrite rod, L, a tuning capacitor, C, and a secondary winding for 
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impedance transformation from antenna to transistor input. An input 
equivalent circuit of Figure 2a is shown in Figure 2b where e is the 
voltage induced by the signal field, r is the series resistance of the 
tuned transformer, and is the transistor input resistance. Another 
equivalent circuit is shown in Figure 2c which is in the form of Figure 




OSCILLATOR INPUT 



o) TRANSISTOR MIXER STAGE 



■or 



b) EQUIVALENT INPUT CIRCUIT 



'a* U. "^rNOOCEO 




cj REVISED EQUIVALENT INPUT CIRCUIT 




d) EQUIVALENT TUBE INPUT CIRCUIT 

Fig. 2 — Mixer stage utilizing loop antenna. 

lb. The voltage V^^, corresponding to a 20-decibel signal-to-noise ratio 
is given by 



10 \/4kT~R^KfjF' 



0.3 



(9) 
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which is derived in Appendix II. In terms of the induced voltage, 
Equation (9) may be written 



\Q^JAkTR,^f^F 

^inducod = • (10) 

0.3 Q, 

The signal input is usually described in terms of field strength so that 
making use of Equation (10), 



Field Strength ^ — = ■ (11) 

h 0.3 Qo h 

for a 20-decibel signal-to-noise ratio, where h is the effective height 
of the antenna, 

A comparison of the noise performance of a transistor mixer stage 
and a tube converter stage may be made by the use of Equation (11). 
An equivalent input circuit for a tube stage is shown in Figure 2d 
where the primary of the antenna transformer of Figure 2c is con- 
nected directly to the signal grid. The value of the tube noise factor 
to be used in Equation (11), depending on the input circuit pai^ameters 
and the type of converter tube, may be determined from Equation (7). 
For typical converter tubes, R,^ is of the order of 150.000 to 300,000 
ohms and the corresponding is of the order of 1,5 to 6 decibels. 

Examples of Noise-Perforniance Calculations for Loop Antennas 

Noise-performance calculations will be made for two battery-oper- 
ated receivers employing feri'ite cored loop antennas, one receiver 
having a transistor mixer first stage and the other a pentagrid con- 
verter tube first stage. The following information applies for both 
receivers. The antenna core consists of a ferrite rod 0.25 inch in 
diameter by 7 inches long. The variable tuning capacitance is 8 to 172 
micromicrofarads. The effective height of the antenna is about 0.004 
meter at 1,000 kilocycles.* The noise bandwidth of the receivers is 
3 kilocycles. 

The receiver employing the transistor mixer has an unloaded 
antenna circuit Q of 200. Substituting the pertinent values in Equa- 
tion (11), 



* This value of effective height applies for a particular core material 
and inductance. It should be noted that for a given core the effective 
height increases with increasing inductance, but for a constant Q the ratio 
y^o/Qoh appearing in Equation (11) is unchanged. 
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10\/1.6 X X 612 X 10^ X 3 X 10-'^ F 

Field Strength = ~ = 

0.3 X 200 X 4 X 10-3 

226v"^ microvolts per meter 

for 20 decibels signal-to-noise ratio at 1,000 kilocycles. For a transistor 
converter having a noise factor of 6 decibels, the field strength for a 
20-decibel signal-to-noise ratio is 452 microvolts per meter. 

The second receiver employs a 1R5, a popular converter tube in 
battery portables, having an equivalent noise resistance of about 
200,000 ohms. The antenna unloaded Q is usually limited to about 
100 by other circuit considei^ations. (Signal feedback, loss of high- 
frequency audio response, and loss of sensitivity due to tracking errors 
are all aggravated by increased Q.) Substituting the required terms 
in Equation (11), 

lOVl.O X 10-20 X 306 X 10^ X 3 X 10^ F 

Field Strength = ■ = 

0.3 X 100 X 4 X 10-3 

^20\/F microvolts per meter 

for a 20-decibel signal-to-noise ratio at 1,000 kilocycles. The converter 
noise factor, as obtained from Equation (7) is 2.2 decibels at 1,000 
kilocycles, resulting in a field strength of 412 microvolts per meter for 
a 20-decibel signal-to-noise ratio. 

Insertion Loss and Image Rejection 

Equivalent input and output circuits for the r-f stage of Figure la 
are shown in Figure 3. These circuits are convenient for the consid- 
eration of insertion loss and image rejection. In Figure 3a, the voltage 
source of Figure lb is replaced by an equivalent current source, and 
the transistor input resistance is referred to the transformer primary 
and designated Ri. The transistor input power is 

i\ = r~R^2 ^ (12) 

which maximizes at i^Ri/4 when R^ ~ R^. A curve of input power 
relative to maximum transistor input power (hereafter referred to 
as insertion loss) for varying Ri/R^ is shown in Figure 4. 

The operating Q of the antenna transformer, designated Qj, is 
related to the unloaded Q by the relation 
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1 



Qoi, 



(13) 



Ro 



1 + 



72i 



and the image rejection ratio, /, is very nearly 



/ / fo 

/-Qi ■ 

\ fo f 



(14) 



where fo is the frequency of the desired signal and / is the image 



frequency. Equation (12) is valid^ for values of greater than 10. 
Maximum image rejection is obtained when the transformer is unloaded 
(Qj — Qqj). The loss of image rejection x^elative to this maximum is 
also plotted in Figure 4. 

A significant improvement in image rejection at a small sacrifice 
in insertion loss results from operation at values of Ri/Ro greater than 
unity. For instance, a change from R^/Rq = 1 to R^/Rq = 2.5 results 
in an improvement of 3 decibels in image rejection at the expense of 
less than 1 decibel in insertion loss. 

An output equivalent circuit which consists of a generator coupled 
to a load by means of a tuned transformer is shown in Figure 3b. The 
elements R^i, Rj^, and R^ represent the output resistance of the r-f 



^ F. E. Terinan, Radio Engineers Handbook, McGraw-Hill Book Com- 
pany, Inc., New York, 1943, p. 144. 




(a) EQUIVALENT INPUT CIRCUIT 




(b) EQUIVALENT INTERSTAGE CIRCUIT 



Fig. 3 — Equivalent circuits of Figure 1 r-f stage. 
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stage, the intei^stage transformer tuned resistance, and the load re- 
sistance, respectively. Maximum operating Q (maximum image rejec- 
tion) for a given insertion loss obtains for the condition Rf^ — R^. 
For this condition the power into the load resistance is 



S 5 
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Fig. 4 — Variation of insertion loss and image rejection with the ratio ■ 
Ri/Ro for the antenna circuit of Figure 1. 



R,, \2 
R^ ( 2 — - Hh 1 

and the operating Q of this transformer is 

1 



(15) 



1+2- 



R. 



R. 



(16) 



where Q02 is the unloaded Q. 
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The insertion loss, determined from Equation (15), is shown in 
Figure 5 as a function of Rj^/Rj^. The image rejection relative to an 
unloaded tuned transformer, also shown in Figure 5, is obtained from 
Equations (14) and (16) for varying Rq/Rj^, 

Example of R-F Stage Design 
The design of an r-f stage which provides near optimum noise 















































NSERTM 


DN LOS 
























































































IMAGE 


REJECT 


ON LOS 


S 







































1.0 1.5 2.0 2.5 3.0 J5 4.0 



12 > 



Fig. 5 — Variation of insertion loss and image rejection with the ratio 
Rb/Ra for the interstage circuit of Figure 1. 



performance and 66 decibels image rejection in a receiver employing 
a 455-kilocycle i-f will be considered. The antenna and interstage 
transformers which were used each have an unloaded Q of 70. 

Near optimum noise performance (determined experimentally) is 
obtained in the range of R^/Ro from 0.4 to 2.5; for this example let 
R^/R^-^l.b. From Figure 4, the image rejection is 4.3 decibels below 
maximum. The maximum image rejection as calculated from Equation 
(14) would be 39.7 decibels at 1,000 kilocycles for = 70. Therefore, 



www.americanradiohistorv.com 



TRANSISTOR RECEIVERS 



155 



the antenna transformer provides 35.4 decibels of image rejection and 
an insertion loss of 0.2 decibel over that obtaining for the maximum 
input power condition. The remainder of the required image rejection, 
30.6 decibels, is to be obtained from the interstage transformer. This 
allows for a tolerable loss of image rejection of 9.1 decibels, relative 
to an unloaded transformer. As determined from Figure 5, the ratio 
Rji/R^ is 0.92 and the consequent insertion loss is 3.8 decibels. 

The r-f transistor has matched terminating resistances of 250 
ohms input and 10,000 ohms output and 30 decibels matched gain; 
the mixer input resistance is 500 ohms. The capacitance required to 
tune the antenna transformer is 100 micromicrof arads ; the capacitance 
required to tune the interstage transformer is 500 micromicrofarads. 
From the above information the turns ratios of both transformers may 
be determined. The tuned resistance of the antenna transformer at 
1,000 kilocycles is 

1 

/e^r=Q^Xcj. = 70 = 112,000 ohms. 

277 10« X 10-10 

For the turns ratio, n, of the antenna transformer, 



1.5 = = , or n — 



Ro nm, V 1.5 R, 

Inserting the vahies of R^ and R^^ in this expression, 



250 1 



1.5 X 112,000 26 

assuming unity coefficient of coupling. The turns ratios of the inter- 
stage transformer may be calculated in a similar manner. The tuned 
resistance at 1,000 kilocycles is given by 



Rj^ =z Qo^> X, = 70 = 22,400 ohms. 

27f 10'> X 5 X 10-10 

The turns ratio from the collector to the top of the tuned circuit is 
given by 



Mc / 0.92 X 10,000 1 

Mrn y 22,400 1.56 
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and the collectoi'-to-secondary turns ratio is given 



Ma 



Ms 




10,000 



500 



0.224 



1 



The r-f stage gain from the base of the r-f transistor to the base of 
the mixer will be (30 — 3.8) or 26.2 decibels. The circuit with nor- 
malized turns ratios is shown in Figure 6. 

In the above design procedure, consideration was limited to opera- 
tion at 1,000 kilocycles. This procedure is satisfactory only if perform- 
ance at the extrem.es of the band is not severely compromised. The 
relative change in image rejection obtained with two tuned circuits 
which are either unloaded or loaded by a constant amount (Rq/Ri and 
Rb/Ra constant) is shown in Figure 7a, Curve A. If the transistor 
input and output resistances are constant over the band, the ratios 
Ro/Rx and Rb/Ra ^vill not be constant. Both the insertion loss and 



image rejection will differ, at the extremes of the band, from the 
mid-band value. The insertion loss will increase and the change in 
image rejection decrease, with respect to the constant-load conditions, 
as the frequency increases. The amount of change will depend on the 
design values of RJRi and Rb/Ra and may be obtained from Figures 
4 and 5 with the aid of Equations (11) to (15). 

The calculated variation in image rejection and insertion loss over 
the bi'oadcast band for the circuit of Figure 6 is shown in Figure 7, 
Curve B, assuming constant transistor input and output resistances. 
Note in Figure 7b that while about ±2 decibels change in sensitivity 
results, the ratio of Rq/Ri varies from 0.54 to 1.6 of the mid-band 
value, and is thus maintained within the previously defined limits for 
good noise performance. For a practical case with many present day 
transistors, the input and output resistances may decrease almost as 
rapidly as R^ and Rb- Curves A of Figure 7 will then be a closer 
approximation to the actual performance. 




F'ig. 6 — Example of r-f stage input and interstage design. 
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Minimum Insertion Loss for Prescribed Image Rejection 

While the above design procedure usually yields satisfactory results, 
it is interesting to consider another approach. The total insertion loss, 
Z, in the antenna and interstage transformers as obtained from Equa- 
tions (11) and (14) is 



UJ p 

UJ 0 
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Pig. 7 — Image rejection and insertion loss of circuit of Figure 6: (A) con- 
stant loading; (B) constant transistor input and output resistances. 



Rn 



Ra 



+ 1 



R, 



R, 



(17) 



and the total image rejection, obtained by combining Equations (13), 
(14), and (16), is 



1 



1 



X 



R. 



Ri 



+ 1 



2Kj, 
Ra 



f fo 
fo f 



(18) 
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Fig-. 8 — Minimum insertion loss versus image rejection at 1,000 kilocycles 
for r-f stage of Figure 1. 

The values of Rq/Ri and Rb/Ra which give minimum loss for a given 
image rejection may be obtained from Equations (17) and (18) and 
are derived in Appendix III. A plot of minimum insertion loss versus 
image rejection, as obtained from Equation (17), for various values 
of ^^^d intermediate frequency of 455 kilocycles is shown in 
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Fig. 9 — Ro/Ri and Rb/Ra versus image rejection (these values provide 
minimum insertion loss at 1,000 kilocycles for the r-f stage of Figure 1). 
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Figure 8. A change of 2-to-l in the Q„i Product causes a 6-decibel 
change in image rejection for constant insertion loss. If, however, 
the image rejection is to be held constant, which is the usual case, the 
change in insertion loss for a change in Qoi Q02 depends on both the 
value of image rejection and the intermediate frequency. A curve for 
an i-f of 260 kilocycles is plotted for comparison. Note that the 260- 
kilocycle curve indicates a loss of image rejection of 8.3 decibels (at 
1,000 kilocycles) as compared to the 455-kilocycle curve for the same 
QoiQo2- Corresponding curves of R, and Rn/R,u ^I'om which the 
transformer turns ratios may be determined, are shown in Figure 9. 

Referring to the sample design in the previous section for a 455- 
kilocycle i-f and Q^^ = Q^y.^70, Figure 8 indicates an insertion loss 
of 3.7 decibels for a 66-decibel image rejection. The values of Ro/R^ 
and Rd/Ra, fi'om Figure 9, are 0.505 and 1.05, respectively. Note that 
this value of Ro/Ri is within the defined limits for good noise per- 




Fijr. 10 — Example of r-f stage input and interstage design for minimum 

insertion loss. 

formance. The corresponding turns ratios, determined as in the above 
example, are shown on the circuit diagram of Figure 10. For this 
design a negligible decrease in insertion loss (0.2 decibel) is obtained; 
other values of QoiQo2 and image rejection may provide a more striking 
difference. This method produces near optimum noise performance 
only by chance. Therefore, its chief value lies in evaluating the results 
of the design compromises of the previous section. 



Appendix I— Derivation of Signal Power Delivered by a 
Capacitive Antenna 

Referring to the equivalent circuit of Figure lb and using the 
notation of the text, the signal voltage at the primary of the trans- 
former is given by 
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where is the loaded (operating) Q of the transformer. The signal 
voltage at the secondary winding is then Vso — '^^a (Ca/Ct) Qi, and 
the signal power delivered to the transistor is 



n 7, Qi 



Since the operating Q is related to the unloaded Q, 



1 

Ri 



R,fii/n- 1 R^ 

Ql = X = Qni 

+ Ri/n^ oyL R^ + Ri 



The signal power may be expressed as 



R. 



Qni 



I Ri' 



(1) 



Appendix II — Derivation of Signal Required for 20 Decibels 

SiGNAL-TO-NOISE RATIO FOR A LOOP ANTENNA 

Referring to the equivalent circuit of Figure 2c and using the nota- 
tion of the text, the signal voltage at the primary of the transformer is 



Ri/n"- 



Ro + Ri/n- 



V., 



Q 
Qo 



The signal voltage at the secondary winding is then Vg, — n Q/Qo> 
and the signal power delivered to the transistor is 

(Vs,)'- 1 

= = in V,, Q/QJ = X — . 

Ri Ri 

The noise power referred to the transistor input is 

Rj 

L Ro -f Ri J Ri 
where F is the transistor noise factor. 



P,v = 4JcT R^Af 



1 
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For an output signal to noise ratio of 20 decibels and 30 per cent 
modulation, 



or, 



= 100: 



4kTn^ R^ Af^ 



_n-R^ + Ri J Ri 



2 1 

— F 



lO^AkTR, Af^,F 



0.3 



(9) 



Appendix III — Derivation of Minimum Insertion Loss for a 
Specified Image Rejection 

Equation (18) of the text is 



Ro 2R, 
— + 1 ■ ^+1 



/ fo 
fo f 



which may be written as 



R, 



+ !=■ 



/ fo V 

fo f 



2Rji 



R. 



+ 1 



Introducing a constant, 



1_ 

fo f 



I 



(18) 



which depends on the required image rejection, the intermediate fre- 
quency, and the unloaded Q's, Equation (18) becomes 



R. 



R. 



■ + 1=- 



Ri\ 



Ra 



+ 1 
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1 



^1 



+ 1 



4 {R,jPh) 



(17a) 



C-[ 1-1 



Differentiating Equation (17a) with respect to R„/Ri, 



R. 
R, 



R„ 



R„ 
R^ 



4 R^ 



(1-C)^ ^2 



+ 2 



Ro 
R, 



+ 1 



(1 C) + 



R,. 
R, 



Ro 
R, 



R„ y 

R, I 



4^^-^ )(1-C) + (1-C)2 

/ R„ Y R„ 

2 (1--C) + — (1 o- 

\ A, / R, 

Equating this expression to 0 to determine the value of Rg/Ri for 
minimum I, 



(C + 2) + \/C^ + 8C 



H 2 
and the corresponding value of Ru/Ha is 

Rj, C 



I 



A',1 - (C + 2) + VC- + 8C 2 
The minimum loss for a given C is then given by Equation (17a) : 



C' [-C + \/C- + 8C] ^ 

2 [- (C + 2) + VC^ + 8C ] [ 3C - VC^ + 86' ] " 
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OPTIMIZING THE DYNAMIC PARAMETERS 
OF A TRACK-WHILE-SCAN SYSTEM 

By 

Jack Sklansky 

RCA Laboratories, 
Princeton, N. J. 

Suntmarij—Systems such as search radar, in which the input data 
arHves intermittently, frequently require a device for eontbmously esti- 
mating the ''present" value of the input. In radar terminology, this device 
is called a ''track-while-scan system" 

A common type of track-ivhilescan system is characteHzed by two 
dynamic parameters, one parameter correcting position error, and the 
other connecting the velocity error. Such a system is analyzed here, and a 
scheme for optimizing the two dynamic parameters is suggested. Tools are 
deHved for the optimization scheme in the form of charts and foi^ulas 
describing the stability, transient response, noise, and manetiver error a^ 
functions of the dynamic parameters. 

One interesting result is a formula for the mean square response of 
the system to white noise. 

A numerical example illustrates a suggested optimization procedure 
using the derived charts and formulas. 

Introduction 



T 



-N^HE efficient use of a search radar in which one or more targets 
appear on the screen intermittently usually demands a device 
for tracking the targets automatically. Such a device, called a 
"track-while-scan system/' must make an estimate of each target's 
instantaneous position from the sampled-data information provided by 
the radar. In a more general sense, the term *'track-\vhile-scan system" 
may denote any system which estimates the "present" value of a signal 
from the "past" sampled values of the signal, the sampling taking place 
at regular intervals. 

In this paper, a common type of track-while-scan system is analyzed, 
and graphs and formulas for optimizing the system's dynamic param- 
eters (or "smoothing constants") are derived. 

In the particular system to be analyzed, the target motion is re- 
stricted to a plane surface, and the estimated course is a sequence of 
straight-line constant-velocity paths, successive corrections in the com- 
puted position and velocity of the target being made proportional to 
the deviations of the computed from the measured positions at sampling 

1G3 
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instants. 1 It is assumed that the data fed to the track-while-scan 
system (the "input") is in cartesian coordinates; the information is 
ordinarily obtained from a coordinate converter operating on the raw 
polar-coordinate data. An illustrative input wave and the resulting 
form of the output are shown in Figure 1. 

The constants of proportionality, « and used in correcting the 
position and velocity, respectively, of the estimated target course 
completely characterize the performance of this particular track-while- 
scan system. These constants are the dynamic pax^ameters or so-called 
"smoothing constants" of the system. The purpose of this paper is 
to develop graphical and analytical tools for optimizing these smoothing 
constants with respect to the following four aspects of performance: 




I ^ 1 1 \ \ 

0 T 2T 3T 5T 6T 

TIME 

Fig, 1— A track-while-scan input and the resulting output. For illustrative 
purposes the input here is varying much more rapidly within each sampling 
period than would occur in practice. 

stability, transient response, noise, and maneuver error. Much of the 
analysis makes use of "z-transforms," a tool especially suited to linear 
sampled-data systems. One interesting result is a formula for the 
mean-square response of the system to white noise. 

A numerical example serves to illustrate how the charts and 
formulas may be used. 

Preliminary Analysis 

In the system under consideration, it is assumed that the target 
path lies in a plane, or that one is interested only in the projection of 
the target onto some flat surface such as the horizontal plane. It is 



1 A less general analysis of such a system has been described by E A 
Mechler, J. B. Russell, and M. G. Preston in "The Basis for the Optimum 
Aided-Trackmg Time Constant," Jour. Frank, Inst., Vol. 248, p 327 Octo- 
ber, 1949. ' 
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also assumed that the input data is received in the form of cartesian 
coordinates x^(nT) and y^{nT), or is converted to this form by some 
appropriate coordinate converter. The track-while-scan system is as- 
sumed to operate on each coordinate separately, and to combine the 
outputs vectorially. 

Successive corrections to the position and velocity of a stored target 
coordinate Xp{t) (henceforth referred to as the output) ai'e made 
proportional to the differences between the obsei^ved coordinate x^{t) 
(henceforth referred to as the input) and the output Xp(t) at the 
sampling instants. As time progresses from the interval (n--l)T < 
t^nT to the next interval, nT < t ^(n -\- 1)T, the .xyaxis intercept 
of Xj,(t) is raised by the amount a[x^{nT) —Xp(nT)], and the slope 
of Xp(t) is increased by the amount plx^inT) —XpinT)'], a and /B 
are the '^position smoothing constant" and "velocity smoothing con- 
stant," respectively, of the system. 

If neither « nor /3 is zero, a constant-velocity target will eventually 
be followed with zero error, assuming there is no noise. However, 
some noise in the input is unavoidable, and the target w^ill not always 
follow a constant-velocity course. Obviously, the system varies in its 
ability to smooth noise and reduce maneuver error for different values 
of a and ^. Techniques are described for determining values of « and 
p w^hich wull minimize noise response and maneuver error while main- 
taining acceptable transient response and stability. 

The track-while-scan system described above is equivalent to a 
sampled-data feedback system in which: (a) the error, defined by 
e^(t) = x^(t) —Xj,(t), is converted to a train of impulses, 

e/(t) ^J^e,(nT)S(t-nT), (1) 

n 

where 8(t) is the Dirac delta function; and (b) this train is applied 
to a single- and double-integrator combination whose impulsive re- 
sponse, g(t), is the sum of a step function of height a and a ramp of 
slope f3/T, The block diagram is shown in Figure 2. The single- and 
double-integrator combination is labeled by its Laplace transform, 
<^{jf(t) ~ a/s p/Ts- ^ G(s) y and the transformation of the contin- 
uous error e^(t) into the impulse train e/ (t) , defined in Equation (1), 
is indicated by the sampling switch. 

Because Xp{t) ordinarily contains discontinuities at the sampling 
instants, it is important to remember that the sampling switch samples 
e^(t) at instants just preceding the discontinuities in Xj,(t). In other 
words the amount of correction which takes place in each new sampling 
period is proportional to the error at the end of the previous sampling 
period, and not to the error at the beginning of the new period. This 
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is because in any physical implementation g{t) will not have a true 
discontinuity at the origin, but will rise with a finite, positive slope. 
Account is taken of this fact in the derivation given below of the system 
transfer function (specifically, at Equation (4)). 

With the aid of the block diagram, the input and output signals 




G(s) 



Fig. 2 — The sampled-data feedback system equivalent to the particular 
track-while-scan system under analysis. 

can be related by Laplace transforms, and the difference-equation 
technique often used in the analysis of sampled-data systems can 
thereby be avoided. 

Using results available in the literature,^ one can relate the output 
and input immediately: 

X,(s)= , (2) 

1 + G*(z) 

where Xp(s) is the Laplace transform of Xp(t), X^'^iz) is the Laplace 
transform of x^(nT)B(t — nT) , with the variable z replacing 

n 

Quantities such as G(s) and G*(z) in Equation (2), and others appear- 
ing elsewhere, are defined similarly. 
Now, since 

Gis)=— + , (3) 

s T.<?- 

it follows that 

CO 

= E + ^-"] - (4) 



2 J. K. Ragazzini and L. A. 5!adeh. "Analysis of Sampled-Data Systems," 
Trans, A.LE.E., Vol. 71, p. 225, November, 1952. This paper includes 
derivations of the basic Laplace and *V.-transform" techniques for the 
analysis of sampled-data systems. 
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The subtraction of a is required here because, as explained above, 
sampling at the discontinuity of g (t) takes place at an instant just 
preceding the discontinuity. 

It can be shown that Equation (4) sums to the following closed 
form: 

(a + 13) z-^ —az - 

G*(z)^ . (5) 

(l-z-"-)- 

Putting Equations (5) and (3) into (2) yields 

aTs + p 



(1- 

A'^(.s) r.- X^;(z) ^ . f6) 

1 - (2 — a~l3)z'-'' + (l~~a)z-- 

This is the desired expression — it completely describes the output .r^, (0 
in terms of the sampled input data x\,(n7')/ This expression will be 
the starting point in the remaining analysis of all four properties of 
interest: stability, transient response, noise, and maneuver error. 

Stability 

From Equation (2) it follows that a necessary and sufficient con- 
dition for stability is that the zeros of 1 + G*(e^^') lie within the left 
half of the plane, or, equivalently, that the zeros of 1 + G*(z) lie 
inside the unit circle of the z plane. From Equations (5) or (6), it 
is seen that the zeros of 1 ^ G* (z) are the roots of thf^ following 
equation : 

z^ — (2 ~ cx — I3)z + (1 ~ a) 0. (7) 

By transforming the interior of the unit circle of the z plane into the 
left half of a new y plane via the transformation y = (z — 1) / {z + 1) y 
one can apply the Routh-Hurwitz conditions for stability to the 
coefficients of Equation (7). The resulting requirements for stability 
are « > 0, /? > 0, and (2a + f3) < 4. An additional stable condition 
exists when ^ = 0, since in that case the zero of the denominator of 
Equation (6) at ;2; = 1 is cancelled by one of the zeros of the numerator. 



* When a = — Equation (6) reduces to a transfer-function descrip- 
tion of the "first-order hold" discussed in Reference (2), Hence the track- 
while-scan system of this paper may be regarded as a generalized first-order 
hold. 
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The resulting necessary and sufficient conditions for the stability of 
the track-while-scan system are a > 0, ^ ^ 0, and (2« + ^) > 4. These 
inequalities determine a "stability triangle" in the a-p plane, for which 
all internal points and all points on the base (/3 — 0) in the intei^-al 
0 < « < 2 correspond to a stable system. This triangle is shown in 
Figure 3. 




0 0.5 1.0 1.5 



Fig. 3— The stability triangle. 

The conditions for underdamped, critically damped, and over- 
damped transient response are found by inspecting the sign of the 
discriminant of Equation (7) if the roots are complex, or the signs 
of the roots themselves if the roots are real.* The resulting conditions 
are 



• If the roots are real, a positive root corresponds to an overdamped 
or critically damped natural mode, while a negative root corresponds to a 
damped natural mode oscillating at one half the sampling frequency. 
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(a + ^) - < 4y5 <— underdamped, 

p^l, {a-\-py- = Ap <=> critically damped, 

^<1, rt^l, {a^p)^>Ap <=:^> overdamped. 

All other values of {a, The transient response contains 

p) inside the stability at least one damped oscillatory 

triangle and on the <=> natural mode with a rate of oscil- 

base {p — 0) in the lation equal to one half the sam- 

interval 0 < a < 2 pling frequency. 

The corresponding regions in the stability triangle are indicated in 
Figure 3. 

It should be noted that the condition (« + jS)^ 4^ alone does not 
describe critical damping. When, under this condition, p > 1, the 
transient will have an oscillatory component, although the envelope 
of the oscillatory component has the appearance of a critically damped 
wave. For the region satisfying both (a -j- p)~> 4^, and a or > 1, 
the transient has a similar characteristic, but with an overdamped 
envelope. 

The relations between the stability triangle and the roots of the 
characteristic equation are described in more detail in Appendix I. 



Transient Response 

The transient response of the system becomes significant when the 
input is switched from one time-shared target to another or when the 
target makes a sharp maneuver. In this section, some curves describing 
the transient response as a function of a and 13 are derived. 

It is convenient to evaluate the transient performance of this system 
by considering the position error in response to a step of velocity. The 
reason for this is that a system which is nominally critically damped 
or overdamped in accordance with the labeling of Figure 3 will still 
yield a position response having an overshoot when the input is either 
a single pulse or a step of position. No position overshoot will appear, 
however, when the input is a step of velocity or acceleration. 

Using Equation (5) and some elementary manipulations of 2;-trans- 
forms,^ one relates the 2:-transforms of e^(t) and x^it) : 

AVW (l"-^-^)2AVU) 
E/(z) = . (8) 

If the input is a unit step of velocity, then x^(t) = 0 or i when t < 0 
or t^O, respectively. Hence Xf,*(z) = Tz~^/ (1 ~ z-'^)-, and 
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E/ (z) ^ ■ ■ . (9) 

1— (2 — a — (1 — 

For the underdamped region of the stability triangle, the above ex- 
pression is the 2;-transform of a time function of the form e-«'sinr^ 
while for the critically damped and overdamped regions the corre- 
sponding time functions have no oscillations. 

Usually, a desirable choice of a and ^ will yield a slightly under- 




Pig. 4 — Contours of constant decrement I'atia, r. 



damped transient response to a step of velocity, i.e., an error response 
whose natural mode is of the form €-«^sint't. An indication of the 
violence of oscillation of this function is given by the ratio of a peak 
value of le""'sinv^l to the peak value immediately preceding. This 
ratio, whose definition is illustrated in Figure 4, will be called the 
decrement ratio, and will be denoted by r. This ratio has an additional 
significance — it is roughly equal to the relative overshoot of e^{t) in 
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response to a step of acceleration, provided that the frequency of 
oscillation of the natural mode is several times smaller than one-half 
the sampling frequency. 

Corresponding to any given pair (a,P) , the value of r can be 
obtained by relating (a,/3) to n/v by means of the denominator of 
Elquation (9). It can be shown that the ratio n/v is, in turn, related 
to r by the formula r = exp (— tt??/?;) . In this way, several contours 
of constant r were computed; they are plotted in Figure 4. These 
contours are useful for both analysis and synthesis, i.e., for either 
evaluating the transient response, or for determining regions in the 
plane corresponding to acceptable transient performance. 



The input coordinate data is subject to corruption by random noise, 
such as (a) range jitter and (b) wander of the target's apparent 
center of reflection. It is desirable to know how the choice of a and p 
affects the degree to which the noise is smoothed or exaggerated by 
the system. 

The following figure of demerit was chosen to describe the ineflficacy 
of the system in smoothing the input noise: 



assuming Xf^(t) is white noise. In words, the noise ratio, p, is defined 
as the ratio of the r-m-s value of the output to the r-m-s value of the 
input. The bandwidth of the input noise is assumed to be finite, but 
much greater than all natural frequencies of the system. 

We shall now derive an expression for p- in terms of a and p. From 
Equation (2) it is seen that o^p(t) may be considered to consist of 
a sequence of responses of a system whose transfer function is 
G(s)/[1 -\~GUz)] to the impulse train 



Since x^^it) is effectively white noise, the successive samples x^inT) 
are completely uncorrelated. Hence the successive responses are un- 
correlated, and the mean square of their sum equals the mean of the 
integrated squares of the individual responses to each impulse. The 



Noise 




x,*(t) ^^xJnT) 8 {t-nT). 



n 
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problem thus reduces to finding the ensemble average of the integrated 
square of 



G(s) 



L 1 + 



(12) 



and dividing the result by the ensemble average of the integrated 
square of the section of x^(t) in the interval nT < t < (n + 1) T, where 
n is chosen arbitrarily. Since the interval is T time units long, the 
integral of x^^it) over the interval is on the average equal to T times 
the ensemble average of x^~(nT), 
Hence 



G(s) 



1 + J 



dt. 



(13) 



Now let 



C(s) 



G(s) 



G*(2) 



[ 



(1-^-0- 



aTs + jS" 



1 



_ 1 - (2-a-p) + {1-a) Z-^\ 



(14) 



Let the expressions in the first and second pairs of brackets be denoted 
by A{s) and B* {z) , respectively. It is convenient to express the time 
function ^-^C{s) in terms of (s) and the coefficients appearing 

in the series expansion of B'' (z) . Thus, 



(15) 



where 



and 



c(t) ^ C(s)^ 
a(t) ^ A(s), 
T.b„z—^B*(z). 



Fortunately, a(t) is bounded in the time domain, thereby affording a 

/• 00 
c^(t) dt. The function a(0 is 
0 

the sum of the doublet and the triangular waveforms shown in Figure 
5. Stated in analytical terms, 
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a+ — t, for Q<t^T, 
T 

fovT<t^2T. 

T 



Let c^{t) for N 0, 1, 2 • • • be defined by 

^ ^. ri for^^O, 

where 1(0 = unit step function = ^ ^ q. 

Then, by Equations (15), (16), and (17), 



(16) 



(17) 









0 T 


2T 


1 


I 
1 
1 




0 T 2T 



Pig, 5 — The two waveforms summing to a(t). 



(18) 



where tx = t — NT. Applying Equation (16) to Equation (18) for the 
interval 0 < t-x < we obtain 

cAh)^h^-i(^-- + l^-^^^'^ +6^v(^a+-^^^-y (19) 



0<t^^T. 



Then, 
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Maneuver Error 

Another important factor governing the choice of the smoothing 
constants is the error incurred when the target follows a curved course. 
This ^*maneuver error" is an additional indication of the poorness of 
performance of the system. 

For convenience of analysis, it is assumed that the target is follow- 
ing a circular arc in the x-y plane at a constant angular velocity w. 
The maneuver error, c\t), will be resolved into and ey{t) . 

The results derived in this section describe the steady-state values 
of the maneuver error. For this reason, the symbols e(0, ^^(0 and 
Cyit) in this section denote steady-state values. 

From the geometry of e^{t) , Cyit), and e(f), it follows that 



(27) 



and that the a: and y coordinates of the true target position, namely 
.r„(0 and 2/^(0, will both be sinusoidal time functions, separated in 
phase by 90 degrees, since the target is assumed to follow a circular 
course at constant speed. Thus, we write 



and 



(28) 



where R = radius of the path. In the above equations, x^{t) has been 
arbitrarily set at zero, and y,(t) has been arbitrarily assumed to lead 
rather than lag x^{t). (These assumptions do not incur any loss of 
generality.) Setting z=:e^^'^ in Equation (8), one obtains an expres- 
sion for the ratio of e^(nT) to x,(nT) in the steady state, i.e., one 
obtains the ratio e,XnT) / (Re^^'^^'l for positive integral values of n, 
the values of 7i being taken large enough to make the transient com- 
ponent of e^ijiT) negligible. Thus, 



e,r(nr) 



R 



1 + G*(€^^n 



(29) 



Since 7/,(0 leads x,(t) by 90 degrees it follows that e/nT) must lead 
e^(nT) by 90 degrees. Furthermore, the vectors e^(0 and ey{t) are 
perpendicular in space. Hence their vector sum, e{t), is circularly 
polarized: it is constant in magnitude at all sampling instants, and it 
rotates in equal angular increments, AO, from sample to sample, the 
rate of rotation, AO/T, being equal to the angular target velocity, o). 



1 
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The magnitude of e(nT) and its angular displacement from the radius 
vector joining the center of the circular course to the target position 
is given by the magnitude and phase of the bracketed quantity in 
Equation (29). The magnitude of the total error, e(nT), is 



R 



1 -h G'(e^"^'0 



(30) 



where the argument of e(nT) has been omitted, because e(7iT) is con- 
stant for all integral n. Setting z = e^'^^' in Equation (5) and substi- 
tuting the resulting expression for G* (e^coT^ \^ Equation (30) gives 
an exact expression for e: 



, (31) 

V[(2~«) coso^r— (2 — « — ^)]2 4-«2sin2a)T 

where R has been replaced by a/oj-, a being the centripetal acceleration 
of the target. 

Ordinarily the sampling rate is such that o)T « tt, in which case 
i^lnoyT r-. a>T and coso>f --^^ 1 — o,-rV2. Then Equation (31) reduces to 





a- — 2/? + a/? ~ 




1 + 













(32) 



when itiT « tt. 



If, in addition, the second term under the radical is much less than 
unity, then the expression becomes still simpler: 



, (33) 

P 

cv2 — 2^ + aP 

when oyT « tt, o)-T- « 1, 

Equations (31), (32), and (33) are the desired results. 
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It is noted that the error in all three equations is exactly propor- 
tional to the target's acceleration. Furthermore, if the sampling rate 
is high enough for Equation (32) to hold, then the magnitude of the 
error is almost independent of the target velocity and is proportional 
to the square of the sampling period. The phase lag of the error be- 
comes negligible when the second term under the radical of Equation 
(82) is sufficiently small. These properties are in conformity with 
physical considerations. 

From Equation (33) it is seen that if the number of samples per 
revolution of the target in its circular course is large, and if the 



OUTPUT OF THE 
TRACK-WHILE -SCAN 
; SYSTEM 




Fig. 7 — Diagram illustrating mechanism by which the steady-state error 
remains essentially unaltered under a change in the position constant, «. 

second term under the radical of Equation (32) is negligible, then 
the steady-state maneuver error is virtually independent of «. The 
mechanism underlying this property is demonstrated by the situation 
in Figure 7. In the steady state the successive values of the error are 
constant in magnitude, and the straight-line portions of the computed 
course must therefore take on equal angular displacements in progress- 
ing from one sampling period to the next. Each of these angular 
displacements is equal to A^, the angular displacement per sampling 
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period of the radius vector of the target from the center of curvature 
of the course. 

The first three periods of the figure illustrate a steady-state condi- 
tion for the case ^ = 0.5; the angle between and A^B^ equals 
the angle between A^B^ and A-^B^. At the sampling instant numbered 
n + 2, the value of a is suddenlj^ raised to 1. This brings the error 
in the computed position down to zero at this instant. The angular 
orientation of C2D2 is equal to that which A2B2 Avould have had if « 
had not been changed, because the error E at the end of the preceding 
sampling period equals the steady-state value attained when a was 0.5. 
Therefore, D^C^^ is nearly equal to E/2. This jneans that the velocity 
correction (practically all of which for this figure is y-oriented) for 
the interval beginning with the (n -f 3)^^ instant is only half that 
required to maintain a steady-state condition. Therefore the angle of 
C:iD.^ with respect to the target course is about 50 per cent larger than 
that of C2D2. The error D3C4 is thus roughly equal to 3E/4. Hence 
the velocity correction for the following interval is again too small for 
the steady-state condition, and D4C5 is roughly equal to lE/S. In this 
way the succeeding values of error at the end of each interval approach 
E monotonically until the steady-state value of E is attained. (More 
specifically, the error is roughly (1—2-^) E when A: sampling periods 
following the change in a have elapsed.) Thus, under the assumptions 
of high sampling rate and negligible value of the second term under 
the radical of Equation (32), which are implied by the geometry in 
Figure 7, it is seen that a change in a has practically no effect on the 
steady-state value of e. These assumptions often are realized in prac- 
tice, especially when the targets are relatively slow. 

It is instructive to deduce Equation (33) directly from the simpli- 
fied situation in Figure 7. In the steady state, the difiference in angle 
between a straight-line path such as A^B^ and the path A^B^ of the 
preceding sampling period is equal to the angular displacement A9 
of the target's radius vector during a full sampling period. This 
angular displacement is also nearly equal to the correction in slope, 
namely, divided by the arc length per period, RAO, Thus: 

RAO 

RAO'^ Ro>-T^ 

and e~ = . (34) 

Since 7?a>2 = a, this result is equivalent to Equation (33). 
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Procedure for Optimizing the Smoothing Constants 

The selected values of a and p should as nearly as possible result 
in a system which minimizes both noise and maneuver error, and 
provides dynamic stability and good transient response. Usually a 
three-way compromise will need to be made among noise, maneuver 
error, and transient response. 

In the following numerical example, a procedure for choosing « 
and P from the derived results will be illustrated. For special require- 
ments, variations of the procedure should suggest themselves. 

Numerical Example 

Suppose the target moves in a circular path in the x-y plane with 
a velocity of 400 miles per hour and a centripetal acceleration of Ig. 
Suppose in addition that the random noise in the observed value of 
target position has an r-m-s value of 50 feet, that successive samples 
of the noise are completely uncorrected, and that the sampling period 
is 3 seconds. 

Problem: Find values of « and p which will provide dynamic 
stability and good transient response, and will yield optimally small 
noise and maneuver error in the computed position of the tai'get. 

Solution : Let the optimization criterion be in the form of a figure 
of demerit, F, defined as the sum of the square of the maneuver error 
and the mean square of the output noise. Thus, 

F ^ + (35) 

where 7^ is the r-m-s value of the input noise. The constants a and /3 
will now be chosen so as to yield a minimum value of F inside the 
stability triangle. 

First we compute the value of o)T to determine whether Equations 
(32) or (33) are applicable: 

aT (32.2) (3600) (3) 

oyT = = r= 0.1647 radian. (36) 

V (400) (5280) 

Since o>T satisfies the condition 0)7 « tt. Equation (32) should be a 
satisfactory approximation. For convenience it shall be assumed that 
the values of or and p to be chosen will be such that Equation (33) 
will also be acceptable. We then are able to write in terms of p: 



37,300 

■ ft-. (37) 
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In the a-P plane the contours of constant e- coincide with contours of 
constant jG; thus, in Figure 6 they are horizontal straight lines. It is 
therefore evident that the minimum value of F must lie at a point of 
tangency between one of the constant-e- lines and one of the constant-p- 
contours. To determine this point, the values of F for several points 
of tangency are computed, and the point of minimum F determined 
by interpolation. 

A typical computation of F runs as follows: Consider the contour 
of p~ = 3. Its point of tangency with a constant-/? line occurs at 
^ = 1.63. Hence 

37,300 

e2 — . — 14,030 ft", 

(1.63)2 

p2ji.2^ 3 (50)2 =z 7,500 ft2, 

F=e^ + phi^ = 21,530 ft2. 

By repeating a few such computations, a minimum value of F = 20,100 
ft2 is found to occur roughly at (a,/3) = (0.5,1,93), corresponding to a 
noise ratio of p^ = 4. Since both inequalities of Equation (33) are 
satisfied here, the approximation assumed in Equation (37) is justified. 

Unfortunately, the fact that the smoothing constants just deter- 
mined will give optimum performance for the steady-state tracking 
of a circular course does not guarantee good transient performance. 
The latter is important, because the target often will not follow a 
circular course for more than a 90-degree arc. It is still more important 
when the system switches from one target to another on a time-sharing 
basis. Satisfactory transient performance can be obtained by confining 
one's choice of (a,jS) to the region where position error in the response 
to a step of velocity will be critically damped or will be underdamped 
with a decrement ratio no greater than, say, 10 per cent. This is 
indicated by the shaded region in Figure 8. 

When the figure of demerit, F, of the numerical example is mini- 
mized within this shaded region, the optimum smoothing constants 
occur w^hen ^ is a maximum, namely, (a,j3) = (0.95, 1.28), correspond- 
ing to which — 3.3 and F::^ 31,050 ft^. The reason is that the noise 
component of F is much smaller than the maneuver-error component 
everywhere within the shaded region. 

Thus, in situations where the maneuver is the major source of 
error inside the shaded region of Figure 8, the pair (a:,jS) = (0.95, 
1.28) will be an optimum. Such situations often arise in practice. 

If desired, the permissible values of (a,f3) may be confined to the 
contour of critical damping. In that case, if the noise component in 
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F is smaller than the maneuver error along- this contour — as it is in 
the numerical example — the resulting optimum parameters will be 
Ka,p) = (1,1)- This corresponds to /)- — 2.67 and, for the numerical 
example, F = 44,000 ft^. 




Fig. 8 — The shaded region indicates the values of (a, P) for which the 
decrement ratio of the position error in response to a step of velocity lies 
between 0 and 0.1. In the illustrative example, this region is assumed to 
correspond to a satisfactory transient response. 



Nonmaneuvering Targets 

Although the preceding example involved a maneuvering target, 
the optimization procedure illustrated in it could be applied to a non- 
maneuvering target, i.e., a target flying a straight-line course at 
constant speed. It is intuitively clear that the optimum values of a 
and p for such a target should be significantly different than for a 
maneuvering target, since greater smoothing will usually be desirable. 

In choosing the smoothing constants for this case, it is important 
to take account of a small unintentional maneuver error or "drift" 
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in the target motion. When this is done, the optimum value of a will 
be relatively small, usually below 0.5. 



Conclusions 

With the aid of the graphs in Figures 4 and 6 and the formulas 
given by Equations (25), (31), (32), and (33) , optimizing the smooth- 
ing constants becomes a reasonably straightforward matter. 

The results of this paper are restricted to two-dimensional tracking 
with no more than two independent dynamic parameters, although an 
extension to the three-dimensional two-parameter case may not be too 
difficult, especially if the target can be assumed to move two-dimension- 
ally over several sampling periods. 
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Appendix I 

The regions of the stability triangle labeled in Figure 9 are de- 
scribed below in terms of the roots in the ;<;-p]ane of the characteristic 
equation (Equation (7)) and the frequency of oscillation of the 
natural modes. 

These results were obtained from the following general expressions 
for the roots of Equation (7) : 



± V (2-«-^) ^ — 4 (1 - ^0 



2-a-/3 ± V(« + yS)--4^ 



(38) 



(39) 



Region I: Complex roots with positive real parts, corresponding 
to an underdamped mode oscillating at a frequency below 
1/(4T). 

Region II: Complex roots with negative real parts, corresponding 
to an underdamped mode whose frequency of oscillation, /, lies 
in the interval 1/(4T) </< 1/(27). 

Region III: Two negative real roots, corresponding to two damped 
modes both oscillating at the frequency /=1/(2T), 
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Region IV: Two real roots of opposite sign, the negative root 
corresponding to a damped mode oscillating at the frequency 
/ = 1/(22"), and the positive root corresponding to an over- 
damped mode. 

Region V: Two unequal, positive real roots, corresponding to two 
overdamped modes. 

Arc AB: Two equal, positive real roots, corresponding to critical 
damping. 




Fig. 9— The stability triangle broken down into regions corresponding to 
the forms of the roots of the characteristic equation. 



Line segment BC: Two conjugate, pure imaginary roots, corre- 
sponding to an underdamped mode oscillating at / = 1/(47). 

Appendix II 

Proof of Equation (22), an important relation used in the deriva- 
tion of the formula for the noise ratio, is given below. 

It was stated that if B*{z) is a rational function of z, and if all 
of its poles lie inside the unit circle, then 
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Residue n 

E_ inside 
f^.v-m^.v — unit 

n=0 • 1 

Circle 



(22) 



Proof: B^z)B*{z-^) = 11 ^j^^ 



CO 00 



Let 7n = i — j. Then 



i=0 m = — CO 



(40) 



= L 2-", 

WI = — 00 



where r^= ?>i_m^. 

i=0 

By Laurent's theorem, 



^m== ^ B*(z)B* 



(41) 



(42) 



where the path of integration is a closed curve separating the poles 
of B*(z) and of B*(z"'^). From this integral, Equation (22) follows 
immediately. 
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CONCENTRIC-SHEAR-MODE 455-KILOCYCLE 
ELECTROMECHANICAL FILTER 



By 

R. W. George 

RCA Laboratories, 
Princeton, N. J. 

Summary — This paper descHbes an experimental 455-Jdlocycle electro- 
mechanical filter of simplified design. The filter consists of four magneto- 
strictive ferrite disk resonators operating in the concentric shear mode. 
The entire filter is only .1 inch long and % inch in diameter. Broad-band 
electrical terminations and low insertion loss make this filter ideal for use 
in transistoHzcd equipment. 

Introduction 

TT^LECTROMECHANICAL filters are known to be particularly 
suitable for use in single-sideband communications equipment. 

^ These filters meet the precision performance requirements and 
are smaller, lighter, lower in cost and in some cases more rugged than 
the equivalent L-C or crystal-type filters. There is a much larger 
potential use of electromechanical filters in relatively low-cost equip- 
ment such as automobile and broadcast receivers. The performance 
requirements here are much less critical, but the cost is of paramount 
importance. The work described here w^as carried out with these factors 
in mind. The result is a small filter which is well-adapted to use with 
transistors. While it can be made much more cheaply than currently 
available electromechanical filters, further cost reduction is necessary 
before the filter can be used in low-cost broadcast receivers. 

The filter is an assembly of magnetostrictive ferrite disks resonant 
in the concentric shear mode^ and quarter-wave torsional couplers. A 
complete filter assembly removed from the case is shown in Figure 1 
and the parts are shown in more detail in Figure 2. Separate couplers 
are cemented between each disk. Magnetic bias is obtained from re- 
sidual circular magnetism obtained by passing a current through a 
hole in the center of the disk. Excitation is by application of the r-f 
field radially to the side of the disk. This makes it possible to measure 
the resonant frequency of each disk as well as to use the end disks 



^ W. van B. Roberts, "Some Applications of Permanently Magnetized 
Ferx'ite Magnetostrictive Resonators," RCA Review^ Vol. 14, p. 3, March, 
1953. 
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Fig. 1 — Filter assembly removed from case. 



of the filter for efficient electromechanical transducers. Uniform cir- 
cular bias and exciting field give negligible excitation of spurious modes 
of vibration. 

The transducers and the mechanical system have very low loss so 
that the important power losses occur in the coils which are coupled 
to the transducers. The coils are enclosed in low-loss ferrite cores, 
Figure 3. These ferrite cores minimize eddy-current loss in the shield 
can and concentrate the r-f field in the transducer disk to obtain maxi- 
mum electromechanical coupling to the enei^gized portion of the trans- 
ducer. Broad-band electrical termination is possible when the electro- 
mechanical coupling coefficient is relatively large, and the magneto- 
strictively energized portion of the end resonator, also a transducer in 
this case, is relatively small.^ It will be apparent that electrical damping 
should be in the external terminal load resistance as much as possible 
so as to minimize electrical transmission loss. 

An axial pivot is fixed to each end of the disk-and-coupler assembly. 
These pivots fit into cup-shaped bearings in the two coil cores. The 
cores are in turn supported by ring spacers and the shield can. The 
pivots should not only space the core from the disk but should have 
low electrical loss as well. Dural pivots were found to lower the Q of 




Fig. 2 — Filter subassemblies and parts. 



2 R. W. George, ''Electromechanical Filter for 100 KG Garrier and 
Sideband Selection," Proc. IM.E., Vol. 44, p. 14, January, 1956. 
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the coil appreciably. Cai'e must be taken to ensure that supporting and 
spacing parts do not contribute appreciable mechanical damping to 
the end resonator disks. Difficulty was encountered from this source 
when pressure was put on rounded lucite pivots; this resulted in a 
5 to 10 decibel increase in the insertion loss. Mechanical damping has 
been satisfactorily reduced by making the pivot bearings hard and 
smooth with a coating of A-2 cementl" and using pointed lucite pivots. 
Spacing between the ferrite core and the disk should generally be as 
small as practical, on the order of 0.01 to 0.02 inch. 




FERRITE DISK AND COUPLER VARIATION 

Fig. 3 — -Mechanical filter details. 



Design 

It is possible to calculate the filter dimensions to give the required 
bandwidth on the basis of the kinetic energy of the resonators and the 
couplers. However, it is more practical to use empirical methods. 

An attempt was made to use coupling coefficients of the Campbell 



t Product of the Armstrong Products Co., Warsaw, Ind. 
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type.^ This was tried by using identical couplers with end disks % the 
thickness of the interior disks. Usable filters were obtained; however, 
the attenuation on the low-frequency side of the passband w^as sloppy 
and included unexpected spurious responses. Filters made using iden- 
tical disks with somewhat larger end couplers eliminated these difficul- 
ties. The spurious responses of these filters were generally at least 
60 decibels down over the frequency range of 160 to 1,600 kilocycles. 
One unexplained response about 50 decibels down has been found in 
some filters around 700 kilocycles. 

In the narrower-band filters, the couplers are more delicate and the 
assembly with cement requires more care. For a bandwidth of some- 
what less than 2 kilocycles, it might be better to use thicker resonator 
disks to increase the required mass of the couplers. The thickness 
should probably not exceed about 0.075 inch which is approximately 
a quarter wavelength in the torsional mode. 

The input and output damping or load resistance values can be 
varied over a fairly wide range wath various values of coupling-coil 
inductance and high-Q coils. Low transmission loss is also obtained 
by having high Q, L-C circuits. 

Magnetostrigtive Ferrite Resonator Disks 

The disk resonators were made of nickel-iron ferrite containing 
vei'y little cobalt impurity.'* The resonant frequency of the disk is 
principally a function of diameter. The hole in the center which is 
necessary in order to magnetize the disk contributes a minor modifica- 
tion of the resonant frequency. Nonuniformity in the thickness of a 
disk also results in variations of resonant frequency. Experience with 
small lots of 25 to 50 disks pressed by hand indicates that the 455- 
kilocycle disks might be expected to have resonant frequency variations 
of less than ±500 cycles. In the filters described, it was found that 
the resonant frequency of the completed filter was some 1,100 to 1,300 
cycles lower than the resonant frequency of the disks themselves be- 
cause of the method of assembly of couplers and disks. 

The disks used in the experimental filter had a diameter of about 
0,507 to 0.509 inch and a hole diameter of 0.077 inch. Nominal thick- 
ness was about 0.062 inch. They w^ere selected or tuned to be matched 
in frequency within about ±50 cycles for each filter. It is difficult to 



3 W. van B. Roberts and L. L. Burns, "Mechanical Filters for Radio 
Frequencies," RCA Review, Vol. 10, p. 348, September, 1949. 

^ R. L. Harvey, "Ferrites and Their Properties at Radio Frequencies," 
Proceedings of the National Electronics Conference , Vol. 9, p. 287, February, 
1954. 
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lower the resonant frequency of a disk in the concentric shear mode 
because this requires a removal of mateinal to make the disk thinner 
in the region of the nodal ring zone which has a diameter of about 
0.7 times the disk diameter. If the disks are to be tuned, it is more 
practical to make them a little low in frequency and then raise the 
frequency by grinding the outside edge or corners. 

Physical faults such as voids and cracks not only reduce the electro- 
mechanical coupling coefficient and residual magnetic bias, but also 
result in spurious resonant frequency responses; the main i^esonant 
frequency is usually altered considerably and is accompanied by lower 
amplitude at resonance. 

The temperature-frequency coefficient of these filters is of the order 
of 40 parts per million per degree centigrade. This coefficient can be 
decreased by improving the characteristics of the ferrite. 

The assembled mechanical portion of the filter was re-magnetized 
by the standard current to insure against loss of bias due to shock or 
exposure to strong magnetic fields while handling and assembling the 
disks. It may also be desirable to age the filter at an elevated tempera- 
ture to stabilize the residual bias if operating temperatures are ex- 
pected to be higher than usual. 

Data on Experimental Filters 

The experimental filters were first tested in a completely closed 
shield with leads at each end. Stray input-output coupling was very 
small, being more than 60 decibels down. This coupling is minimized 
by the use of small coupling coils and cores which have an outside 
diameter considerably less than the ferrite disk diameter. The stray 
coupling obtained with the shield can shown in Figure 1, having all 
leads at the open end of the can, was somewhat more than 40 decibels 
down. A material improvement results from closing the open end of 
the can. The ferrite coil cores may be increased in size somewhat to 
increase the coupling which might be required for very-wide-band 
filters, or to make more space for more efficient or higher impedance 
coils. 

The circuit shown in Figure 4 was used in testing these filters. 
Some variation in coil-core spacing to the filter resulted in unequal 
values of R^ and R^ and and Co. Tuning capacitance was on the 
order of 500 micromicrofarads when using 75-turn coils. Tests were 
made with an output voltage, £7o, of about 0.1 volt; however, the filters 
operate satisfactorily with output voltages of up to 10 volts or more. 
Response curves for 4 filters are shown in Figures 5 to 9. The shape 
factor is on the order of 3 to 3.5. 
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SHIELD CAN 

Pig. 4 — Filter test circuit. 




440 450 460 470 

FREQUENCY (Kc.) 



Fig. 5— Filter No. 8. 
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440 450 460 4 70 

FREQUENCY (Kc.) 

Fig. Filter No. 9. 
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FREQUENCY (Kc.) 

Fig. 7— Filter No. 10. 



Electrical insertion loss is estimated to be less than 6 decibels. 
This was based on the power input to the filter at the terminals where 
El is measured, and the output power in i?2- varies considerably 
with frequency in the pass-band with the result that this method of 
measurement is not very satisfactory. The coils were layer or random 
wound with about 75 turns of No. 36 enamel wire. Lower insertion loss 
could be achieved by use of higher-Q coils. 

Mechanically perfect joints between resonators and couplers are 
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432 440 450 460 47X) 4»0 

FREQUENCY, (Kc) 

Fig:. 9— Filter No. 12. 

desired. Satisfactory results were obtained with A-2 cement. 

The use of quarter-wave couplers permits the use of the largest 
possible diameter which will give the desired coupling coefficient. The 
couplers should be made of material having low mechanical losses. This 
material should also have a reasonably low frequency-temperature 
coefficient if the ferrite resonator disks have a very low frequency- 
temperature coefficient. Bandwidths obtained with various parts dimen- 
sions are shown in Table I and Figure 3. 

Pivots or other supporting-spacing means which are in the flux 
path of the coils should have low electrical loss as well as good mechani- 
cal characteristics. Pivots should be either loose or cemented in the 
filter element in such a way as to minimize mechanical losses. Cemented 



Table I — Data On Experimental Filters 
Coxip]ers 



Filter 
Number 


Bandwidth 
(ci>s) 


Interior 

d di 
(inches) (inches) 


End 

d 

(inches) 


di 

(inches) 


(ohms) 


(ohms) 


8 


8260 


0.095 


0.062 


0.100 


0.062 


9000 


18000 


9 


9500 


0.095 


0.058 


0.100 


0.058 


18000 


18000 


10 


5340 


0.090 


0.058 


0.095 


0.058 


18000 


9000 


11 


4000 


0.086 


0.066 


0.090 


0.066 


50000 


4000* 


12 


17000 


0.107 


0.058 


0.115 


0.058 


40000 


18000 



* Across a tap on the output coil. Or this tap can be connected directly to a transistor circuit. 
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pivots may appreciably detune the end resonator and therefore should 
have small mass. Another mounting arrangement which was tried was 
a loose-fitting bakelite shaft through the filter. Small washers at each 
end spaced the filter with respect to the coil cores. This method of 
mounting is satisfactory if very little or no pressure is put on the 
spacers so that the filter element is positioned with substantially no 
binding. 

A variation in disk and coupler shapes is shown at the bottom of 
Figure 3. A recess on each side of the disk centers the coupler and 
pivot-bearing assembly. The couplers are simply cut-off pieces of 
tubing with a permissible length variation of perhaps 2 or 3 per cent. 
The outer diameter is kept constant to fit the recessed hole in the 
ferrite disk, while the inner diameter is chosen to give the desired 
coupling coefficient. The couplers may be made of a variety of materials 
including possibly some ceramics. The dimensions given here apply 
only to the use of dural. 
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VAPOR PRESSURE DATA FOR THE 
MORE COMMON ELEMENTS 



By 



Richard E. Honig 



RCA I laboratories, 
Princeton, N. J. 



Summary — The most recent data available coyicerning vapor pressures, 
melting and boiling points, and heats of sublimation have been selected, 
tabulated, and plotted for 57 elements, many of which are of special interest 
to ivorkers in the fields of electronics and high-vacuum technique. It has 
been found convenient to present vapor pressure data graphically as plots 
of log p (mm Hg) versus log T C^K), and also to tabulate the absolute 
temperatures for fixed pressures. This collection contains dctta published 
or available before March 1, 1957. 



N 1948, R. R. Law^ published vapor pressure data for some 30 



substances, presenting all information graphically on a single sheet 



as curves of log versus log T. La\v*s collection is based entirely 
on prewar data. A review article published by Brewei"^ in 1950 contains 
a substantial amount of information obtained during the w^ar years. 

In the intervening years, a large volume of new information has 
come into existence, and the need for an up-to-date collection of vapor 
pressure data has become more and more pronounced. This review, 
begun in December 1956, w^as initially based largely on the 1955 edition 
of ''Metallurgical Thermochemistry" by Kubaschewski and Evans^ 
as well as a comprehensive literature search carried out by the author. 
Upon completion of the first version, two unpublished compilations of 
data, by Stull and Sinke"^ and by Hultgren,^ became available which 



1 R. R. Law, Vapor Pressure Data for Various Substances (A Graphi- 
cal Presentation)," Rev. Sci. Instr., Vol. 19, p. 920, December, 1948. 

- L. Brewer, The Chemistry and Metallurgy of Miscellaneous Materials : 
Thermodynamics, McGraw-HiH, New York, N. Y., 1950, p. 13. 

2 0. Kubaschewski and E. L. Evans, Metallurgical Thermochemistry , 
John Wiley & Sons, Inc., New York, N. Y., 1956, 2nd Ed. 

4D. R. Stull and G. C. Sinke, "Thermodynamic Properties of the 
Elements," Advances in Chemistry, Series 18, American Chemical Socictv, 
Washington, D. C., 1956. 

s R. R. Hultgren, "Selected Values for the Thermodynamic Properties 
of Metals and Alloys," Report of the Minerals Research Laboratory, Uni- 
versity of California, Berkeley, California, 1956, 
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contained critical evaluations of data up to and including 1956. The 
present review is largely based on these three collections. 

Data 

Prei^entation 

Vapor pressure data may be presented very succinctly with the 
help of the equation 

logioP = A r-i + /? \o^,^T ^CT + D (1) 

where 

p — pressure, expressed in this paper in millimeters of 
mercury, 

T— absolute temperature in °K, 

A, B, C, D ^ coefficients characteristic of the element. 

In this system, employed by Kubaschewski and Evans, ^ the coefficients 
are readily tabulated in a minimum of space, but the evaluation of 
pressures for a given set of temperatures is quite cumbersome. The 
converse process — finding temperatures corresponding to given pres- 
sures — can only be done by trial and error and usually requires the 
service of a computing machine. 

A more useful way of presenting vapor pressure data is to tabulate 
pressures corresponding to fixed temperatures, as was done by Stull 
and Sinke,^ or to quote temperatures for given pressures. Both systems 
of tabulation were employed by Hultgren.^ In general, we are inter- 
ested in a fixed pressure range (10~8 < p < 760 mm Hg) which is the 
same for all elements, while the corresponding temperatures may lie 
anywhere between 0 and about 5,000 °K. If curves are to be drawn 
from such tabulations, it is found that the points will be suitably spaced 
if pressure decades are employed. For this reason, temperatures have 
been tabulated at fixed pressures in the present paper. 

Vapor pressure curves have been presented in the past usually as 
plots of log p versus 1/T, largely because the resulting curves are 
straight lines over a limited temperature range and have slopes which 
yield heats of vaporization directly. Unfortunately, this type of tem- 
perature scale compresses the high temperature end so much that it 
is not feasible to present many vapor pressure curves on the same sheet 
over a wide temperature range. From a practical standpoint, the log p 
versus log T plots first introduced by Law^ are by far the most suitable 
and are employed in this paper. 
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Selection and Treatment 

A compilation of this type is based on the data of many workers 
who measured vapor pressures by different techniques under widely 
differing- experimental conditions. Thus it is often very difficult to 
choose between two apparently equivalent sets of results, and at times 
the preference expressed may be subjective and arbitrary. Still, every 
effort has been made to include all data, based on reliable experimental 
measurements, for those elements that are of some practical impor- 
tance. The only exceptions are iridium, palladium, and rhodium for 
which there exist to date only estimates. 

The vapor pressure data available from the three major collections^"^ 
were compared graphically by plotting log p versus log T. Kuba- 
schewski and Evans' data,^ given as the coefficients of Equation (1), 
were first converted into a tabulation of temperatures coi'responding 
to fixed pressures with the help of a computing machine. They were 
then entered on the log p versus log T plots, and smooth curves drawn 
through the points. In the very complete collection by Stull and Sinke^ 
vapor pressures are tabulated as logio P (atm) at fixed temperatures. 
These could be put directly on the log p versus log T plots with the 
help of a special scale. Hultgren's data^ for selected metals, tabulated 
as temperatures at fixed pressures (in atm), could be plotted directly. 

Data based on Kubaschewski and Evans' collection^ and on a liter- 
ature search made by the author had already been plotted when the 
other two compilations became available. The curves were left un- 
changed wherever the pressures agreed with the other two sources to 
within about 20 per cent, which is the accuracy to be expected from 
vapor pressure measurements. Where there were larger discrepancies, 
preference was usually given to the evaluations of Hultgren^ and Stull 
and Sinke.* It was felt that Hultgren's evaluation of a few selected 
metals was likely to be somewhat more thorough than Stull and Sinke's 
collection covering most of the elements, and that either was preferable 
to the extensive but less up-to-date collection by Kubaschewski and 
Evans who treated most of the elements and many compounds. From 
the curves mentioned above, the temperatures corresponding to fixed 
pressures were read off and tabulated. 

For those elements having two or more important gaseous species 
of known concentrations, total vapor pressure curves were obtained 
by adding up the individual curves graphically. These are identified 
by a 2 preceding the chemical symbol. In those cases where the atomic 
species is known to be the predominant contributor to the total vapor 
pressure, contributions from molecular species were neglected and the 
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symbol given without any subscript. For the few elements that consist 
largely of one molecular species, the appropriate subscript has been 
added. 

Results 

Table I contains the desii'ed vapor pressure data, as well as the 
melting and boiling points, literature references, and in many instances 
the temperature range over which the original experimental data were 
obtained. The uncertainty in pressure considered tolerable, ±20 per 
cent, corresponds to a temperature uncertainty of between 1 and 2 
per cent. To simplify matters, temperatures have been tabulated to 
±r between 100° and 1,000''K, to ±5'' between 1,000° and 2,000°K, 
and to ±10° above 2,000°K, corresponding to the smallest increments 
that can be read on the graph. Thus it should be remembered that the 
last figure quoted is not always significant. In the reference column, 
the three collections rather than the original articles have been quoted, 
thereby avoiding the repetition of a very extensive bibliography. 

Figures lA and IB present the vapor pressure data in graphical 
form.* The curves were placed on two separate sheets in such a way 
as to minimize interference. The circled point shown on most curves 
is the melting point. Where the melting point falls outside the pressure 
range of the graph, the letters "s" (solid) or "1" (liquid) have been 
appended to the chemical symbol. The last column of Table I indicates 
on which sheet a given element will be found. 

Table II gives the most recent values available for the heats of 
sublimation, AHoqs, quoted for the individual elemental species at 
298°K. These heats are of considerable practical interest because they 
are a measure of lattice energies. The accuracy of the values quoted 
will lie somewhere between ±0.01 and 10 kcal ^mole, is often difficult 
to assess objectively, and is not available for all entries. For this 
reason, estimate errors have not been tabulated, but indications of 
accuracy may be had from the number of significant figures given. 

Discussion 

Every effort has been made to include in this collection of vapor 
pressure data all information available as of March 1, 1957. 

Aside from the three entries where data are based on estimates 
(iridium, palladium, and rhodium), there are three elements which 



* A limited number of *'wall-size*' reproductions of these charts are 
available; requests for copies should be directed to RCA Laboratories, 
Princeton, N. J. 
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Table II — Heats of Sublimation of Elemental Species at 298°K 



SPECIES 
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56 


7 




99 
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Rb2 


27.6 


4 
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65 
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^ L. Brewer, private communication. 

^ L. Brewer, ''Heats of Sublimation of the Elements,*' University of 
California Radiation Laboratory Report 2854 (revised), February, 1955. 

s J. Drow^art and R, E. Honig, ^*Mass Spectrometric Study of Copper, 
Silver, and Gold," Jour. Chenu Phys., Vol. 25, p. 581, September, 1956. 

^ R. E. Honig and J. Drowart, "Vaporization of Compounds and Alloys 
at High Temperature," Technical Note 1, University of Brussels, Belgium, 
1956. 

'^^ R. E. Honig, "Mass Spectrometric Study of the Molecular Sublima- 
tion of Graphite," Jon7\ Ckem. Phys., Vol. 22, p. 126, January, 1954. 



www.americanradiohistorv.com 



204 



RCA REVIEW 



June 1957 



require discussion: silicon, aluminium and uranium. The vapor pres- 
sure data presented for silicon were obtained by the writer by an 
indirect method." Recent, direct measurements of SiC by Searcy seem 
to have yielded^-^^ a silicon pressure of 1.6 X lO"-*'^ mm Hg at 1,500''K 
(about 10 X lower than the present value) and AHoqs = 113 kcal/mole, 
but to date no further experimental data or details have been received. 
The data for aluminium are based on a study by Brewer and Searcy^^ 
who used BeO and TaC crucibles and observed considerable scattering 
of points. If the TaC results may be discarded, and if the effects of 
BeO-Al interaction are negligible, as suggested by Bi^ewer,^ then 
their data are much more consistent and will yield AH298 ~ 1^-^ 
kcal/mole. The vapor pressures given for uranium in Table I are 
consistent with AHsos^ 117.2 kcal/mole, as quoted in Table II. The 
heat of sublimation preferred by Brewer^ is AH29S = 122.7 kcal/mole. 
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1^ R. E. Honig, "Sublimation Studies of Silicon in the Mass Spec- 
trometer," Jour. Chem, Phys., Vol. 22, p. 1610, September, 1954. These 
silicon vapor pressures are based on a comparison with germanium. The 
values given in Table I of the present report have been adjusted to take 
into account the most recent germanium vapor pressure data (A. W, Searcy 
and R. D. Freeman, ''Measurement of the Molecular Weights of Vapors at 
High Temperature," II **The Vapor Pressure of Germanium and the 
Molecular Weight of Germanium Vapor," Jour. Chem. Phys., Vol. 23, p. 88, 
January, 1955). 

^- A. W. Searcy, private communication. 

^3 L, Brewer and A. W, Searcy, "The Gaseous Species of the Al-AUOa 
System," Jour. Avier. Chem. Soc, Vol. 73, p. 5308, November, 1951. 
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ALLOYING PROPERTIES OF GERMANIUM FREE 
OF EDGE DISLOCATIONS^^ 



By 

C. W. Mueller 

RCA Laboratories, 
Princeton, N. J. 



Sumrnary—Germanhnn that is substantially free of edge dislocations 
has alloying properties ivhich are covsiderably different from those of 
gernianiiini with many edge dislocations. In particular, surface spreading 
in eqidlibrinm processes is difficult to control. Experiments are described 
showing the effect of surface tension and germanium cry st alio graphic 
pla}ie on iridium alloy dot diameter. Alloying with controlled spreading 
may be accomplished by arranging for initial dissolution of the germanium 
by an unsaturated solution. An alloying method is described which consists 
of ivetting an indium dot on germanium at 300°C in hydrogen and then 
using a very rapid heating in nitrogen to the alloy temperature of 550°C. 
The resulting junctions are extremely flat and uniform. 



Introduction 



1 



^HE advances in the techniques of growing- semiconductor single 
crystals have made materials available with a greatly improved 
physical structure. Germanium single crystals that are essen- 
tially free of edge dislocations can now be readily grown. The alloying 
or dissolution properties of indium on germanium that is free from 
edge dislocations is different than on germanium with many edge 
dislocations. 

In the formation of p-n junctions for transistors, it is necessary 
that the area of the junction and its location be accurately controlled. 
The (111) crystal plane is a natural barrier in the crystal that forms 
the desirable planar junctions. ^ If, however, a crystal has no edge 
dislocations, resistance to penetration by the indium from a saturated 
solution is so strong that excessive and uncontrolled horizontal spread- 
ing of the molten dot occurs.^ 

This paper discusses the important factors in alloying on disloca- 



* Presented at The Electrochemical Society Meeting, Cleveland, Ohio, 
October 3, 1956. 

1 0. W. Mueller and N. H. Ditrick, "Uniform Planar Alloy Junctions 
for Germanium Transistors," RCA Review, Vol. XVII, p. 46, March, 1956. 

2 H.^ V. Kettering and J. I. Pankove, '^Crystal Quality as a Factor in 
Controlling the Alloying of Indium onto Germanium Surfaces," presented 
at the I.R.E.-A.I.E.E. conference on Semiconductor Research, Lafayette, 
Indiana, June 26, 1956. 
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tion-free germanium and shows that controlled alloying on dislocation- 
free material may be achieved by a process in which the dissolution 
occurs mainly by an unsaturated indium-germanium solution. 

Physical Model 

The forces concerned with the spreading of a molten dot can best 
be discussed with the help of Figure 1. The temperature cycle of an 
alloy junction is usually adjusted so that a saturated solution of 
germanium in indium occurs, thereby making the process essentially 
hidependent of time. Hence, the total volume of germanium taken in 
;olution is the same for all units. The forces that determine the shape 
of this volume and how these forces can be controlled so that the 
desired shape is always attained will be discussed. 

At the present time an examination of these forces in a funda- 

MOLTEN DOT 




Fig. 1 — Factors influencing spreading of a molten dot. 

mental manner is difficult because of the sketchy state of knowledge 
about the wetting and flow of one metal on another, ^ especially for 
indium and germanium. Since the shape of the solid-liquid interface 
is the result of the combined actions of the several forces, the rates 
of liquid motion in directions "A," '*B,*' and *'C" (Figure 1) are 
important. The rates are not definitely known, but they can be quali- 
tatively discussed and their relative magnitudes illustrated by ex- 
periment. 

The motion of the molten indium can be resolved into horizontal 
and vertical components with respect to the surface of the wafer, which 
is nearly a (111) plane. From examination of junctions we know that 
the fastest horizontal rate is toward the corners of the triangle formed 
by the (111) planes that intersect the surface. Since maximum spread- 
ing must be controlled, the fastest rate must be contained, and this 
is the one considered. 

The *'A" rate (Figure 1) is a horizontal surface rate that depends 
on the surface tension of the molten dot and the wettability of the 

3G. L. J. Bailey and H. C. Watkins, "The Flow of Liquid Metals on 
Solid Metal Surfaces and its Relation to Soldering, Brazing, and Hot-Dip 
Coating,'* Jour. Inst, of Metals, Vol, 80, p. 57, 1952. 
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surface by the dot material. It is influenced by surface conditions, 
temperature, and dot material. 

The '*B" and "C" rates are not influenced by the surface; they are 
mainly determined by the dissolving of the germanium by the indium. 
They are solution or alloying forces and depend upon the temperature 
and the bulk properties of the germanium, i.e., crystallographic direc- 
tion and perfection of the crystal. Tf the molten dot is uniformly 
mixed, the *'B" and "C" rates will be the same when the crystal bulk 
properties are isotropic. Excessive spreading can be prevented by 
decreasing the rates '*A" and ''B" and increasing rate *'C." Experi- 
ments will be described which illustrate how these rates can be greatly 
changed by the choice of experimental conditions. 

In all the experiments, germanium areas free from edge dislocations 
were used. All the germanium wafers were etched for 10 minutes in 
CP4 and then examined microscopically for etch pits. Regions free of 
etch pits associated with edge dislocations^ were then selected. This 
was not difficult in material that had an average count of 100-300 pits 
per square centimeter. 

Indium-zinc balls (0.5 per cent zinc) 0.014 inch in diameter cleaned 
in 1 per cent HCl were used for dot material unless otherwise stated. 
These dots spread slightly less than pure indium dots. Uniform wetting 
was reproducibly attained using hydrogen as the flux at temperatures 
of about 300° C. No liquid fluxes were used. The large undercutting 
of the surface on dislocation-free germanium by the molten dot aids 
in removing any initially unw^et areas. 

Surface Tension 

Several methods are available to reduce rate ''A" which is closely 
connected with surface tension; 

1. The use of a slightly oxidizing atmosphere. 

2. Mechanical constraint by a jig. 

8. Use of a dot material such as lead alloys having a higher surface 
tension. 

In this paper only method 1 will be discussed. Method 2 has the 
disadvantage of mechanical wear at high temperatures and usually 
does not sufficiently limit junction area because the molten liquid can 
undercut the surface. Method 3 was tried and it was found that the 
use of lead inhibits excessive spreading only when lead-rich alloys 
O50 per cent) are used. However, in this case the alloy penetration 
is very shallow. 

S. G. Ellis, "Dislocations in Gennaniuni," Jour, AppL Phys,, Vol. 26, 
p. 1140, September, 1955. 
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The use of a slightly oxidizing atmosphere during alloying is prob- 
ably the simplest way to change rate "A." An experiment that elimi- 
nates surface tension gives an insight into the nature of the bulk 
alloying mechanisms involved in forming a junction. The indium-zinc 
dot is first soldered or wetted down in a pure hydrogen atmosphere 
at 300° C. The alloying step, heating to 550° is then done in a gas 
from which all the oxidizing gases have not been removed — line nitro- 
gen, for example. Traces of oxygen and water vapor attack the indium 
very slightly and form a thin oxide film. This oxide film forms a 




(a) TOP VIEW OF RECRY5TALL1ZED GERMANIUM 
WITH DOT REMOVED- 



6ERWAN!UW 
WAFER 




INDIUM DOT 



^54^— ,5— 1-5-1 ^■^"'^C-f'O^' 

(all dimensions in mils) 

(b) CROSS SECTION OF DOT ABOVE SHOWING 
LARGE UNDERCUT REGIONS- 

Fig. 2 — "Hat brim" type of structure from use of an oxidizing atmosphere 

(not drawn to scale). 



membrane over the dot and restrains the dot in a hemispherical foi'm. 
The diameter of the hemisphere, as determined by a machinist's micro- 
scope or shadowgraph, does not grow during the alloying step. 

The use of an oxidizing atmosphere thus provides good control of 
the surface forces (rate ''A" of Figure 1). However, on close exami- 
nation of the alloy dot, one finds that the germanium surface has been 
undercut, i.e., rate **B" has not been reduced and a **hat brim" type 
of structure has been created as shown in Figure 2. The surface of 
the crystal of Figure 2 is parallel to the (111) plane. Heating and 
cooling was done at a rate of 20° C per minute between 400° and 550° 
C. Typical dimensions of the outer brim, as shown in Figure 2, can 
be large, thus causing high electrical capacitance and large collector- 
to-emitter spacing. 
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Note that the horizontal spread from the original dot is more than 
10 times the vertical penetration. The recrystallized germanium is very 
thin in the outer brim region because as the dot cools, mixing is 
restricted and the amount of germanium available for recrystallization 
is small. One sees that merely restricting the surface rate of spread 
does not by itself restrict the actual junction because of the very slow 
rate at which the indium penetrates the horizontal (111) plane. Con- 
sequently the bulk crystal properties must also be considered as a 
means of coping with the problem of lateral undercutting. 




(100) 



(110) 



(111) 



010" 

Pig. 3 — Cross section of junctions produced on different crystal planes. 

Crystal Axis 

One way of changing directional properties of the alloying process 
is by cutting the germanium crystal so that a plane different from the 
(111) is parallel to the surface. Figui-e 3 shows a series of cross 
section photographs illustrating the results of alloying on the different 
crystal planes when all are processed in the same manner at the same 
time: namely, wetting down of the 0.014-inch spheres in hydrogen 
at 300° C and then alloying to 550"" C in hydrogen with a heating and 
cooling rate of 20° C per minute between 400° and 550° C. Micrometer 
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measurements and the pictures clearly show that there is much less 
spreading on the (100) and (110) planes. Close examination of the 
junctions of Figure 3 as well as the examination of others from which 
the indium was removed, shows that the wetting on the (100) and 

(110) planes was considerably poorer than on the (111) plane. On the 

(111) plane the large undercutting removes unwet areas. This experi- 
ment does not mean that it is impossible to use the (100) or the (110) 
planes of germanium, but that the wetting on these planes is more 
difficult. In other words, a wetting problem is substituted for the 
spreading problem. Consequently, it is desirable to have a method of 
utilizing the Clll) crystal plane. 

Rate ok Dissolution 

If indium in contact with germanium is heated slowly under con- 
ditions approaching equilibrium, then germanium atoms are continu- 
ally leaving, and returning to, the germanium wafer. In this case the 
differences in binding forces acting on the atoms in different crystal- 
lographic locations have a chance to be effective. On the other hand, 
if the solution rate is rapid (as in unsaturated indium) the crystal 
forces are overcome by the more violent dissolving forces and the 
differences of binding energy for variously located atoms have less 
effect. Under a rapid heating, nonsaturated condition, the solution of 
germanium by indium is that of a liquid drop acting at nearly the 
same speed in all directions. Consequently while the liquid is in a 
highly unsaturated condition, the (111) plane does not hinder penetra- 
tion as markedly. Experiments on the solution rate of oriented ger- 
manium wafers in molten indium are consistent with this hypothesis.^' 

At the end of the heating period while the temperature is constant, 
the solution becomes nearly saturated. Hence crystal forces on the 
atoms are comparatively stronger and can level the bottom of the 
junction. 

A very rapid temperature rise promotes germanium solution by 
unsaturated indium. In Table I the large effect of changing the rate 
of temperature rise is indicated. The indium-zinc dots (0.5 per cent 
zinc) on all the units of the table were first soldered in hydrogen on 
dislocation-free germanium at 300° C. The diameter of the dots was 
then measured in two perpendicular directions. An atmosphere of 
line nitrogen was used for the alloying. As shown in the table, the 
units that were heated at 10° C per minute spread an average of 18.7 
mils, while the rapidly heated (1,500" C per minute) units spread only 

5 B. Goldstein, "The Dissolution of Germanium by Molten Indium/' 
RCA Review, Vol. XVIII, p. 213, June, 1957. 



www.americanradiohistorv.com 



ALLOYING PROPERTIES OF GERMANIUM 



211 



TABLE 1 



GROUP 


HEATING RATE 
300°C.TO 550'C. 


MINUTE5 HELD 
AT SSC^G. 


GOOLING RATE 
550*C. TO 300'C. 


FINAL DIAM. 
(MILS) 


INGREASE 
IN DIAM. 

(mils) 


A 


io^c/min. 


3 


iq-g/min. 


35.8 


18.7 


B 


isoo^g/min* 


3 




19.6 


2.3 


C 


I500"g/mIN* 


30 




19.7 


2.4 


D 


125*c/MIN. 


3 


igo'gAiin. 


30-6 


13.4 




NOT CORRECTED FOR HEAT LOSS DUE 


TO THERMOCOUPLE 


LEADS 





2.3 mils. The sidewise alloying spread or undercutting was thus re- 
duced by a factor of about eight. 

In order to determine if the solution reached saturation with the 
very rapid heating rate the top temperature of 550° C was held in 
group C (Table I) for 30 minutes. The data shows that additional 
spreading is negligible. The depth of penetration is about the same. 
Although the rate of rise must be rapid, time at the top temperature 
is not critical. Equilibrium is reached at the top temperature and the 
process is no longer a dynamic process. This is desirable because a 
dynamic process is generally more difficult to control. The decrease 
of temperature from the maximum can be independently controlled to 
give the desired type of recrystallization (see Figure 4). 

In order to secure the fast rate of rise shown, a jig with a large 
heat mass cannot be used because of its thermal inertia. In the experi- 




1- I 

.008" 



Fig. 4 — Cross section of junctions made by using very rapid temperature 
rise: top view, slow cooling; bottom view, rapid cooling. 
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ments described, the germanium pellets were placed on a very thin 
sheet of mica and pushed into a furnace which had been brought to 
the top temperature of 550° C. Laying germanium w^afers on a 
graphite boat and pushing this into a furnace at 550° C produces the 
results shown in group D of the table. Although the spreading is 
somewhat less than for group A, it is not nearly as small as that shown 
by the units of groups B and C with the extremely rapid temperature 
rise. 

A good alloying process should also produce junctions that are 
crystallographically flat. The cross sections of junctions made with 
a very rapid temperature rise as described above are shown in Figure 

4 which shows that the flatness is determined by a crystal plane. These 
units made on the dislocation-free germanium not only show flat 
bottoms but also show sides that are determined by crystal planes. If 

5 junctions are made in a line on a wafer and a cross section of the 
entire group is made, one finds that they are all perfectly formed. 
This degree of uniformity of the junctions made as described above 
is very striking, and shows that, for this type of crystal, a slow rate 
of heating is not as important a parameter as heretofore supposed. 

Conclusions 

By arranging the germanium alloy cycle so that initial dissolution 
by an unsaturated indium solution occurs, the restrictive action of the 
(111) crystal plane may be minimized. This action gives a simple 
method of making junctions without excessive spreading on germanium 
free of edge dislocations. As the final atom layers are removed the 
solution becomes saturated and the restrictive action of the (111) plane 
again becomes very strong leveling the bottom of the junction. A 
two-step process is used. The dot material is wetted to the germanium 
at 300° C in hydrogen. No liquid flux is used. The alloying is then 
done in line nitrogen with an extremely rapid temperature rise to the 
alloying temperature of 550° C. Equilibrium is reached in less than 
three minutes and the temperature decrease can be governed by the 
type of recrystallized region desired. 

The junctions formed by the above process are extremely flat and 
uniform when examined optically at high magnification. 
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THE DISSOLUTION OF GERMANIUM 
BY MOLTEN INDIUM^^ 

By 

B. Goldstein 

KCA Laboratories, 
Princeton, N. J. 

Summary — Experiments have been performed on the dissolution of 
germanium by molten indiam as fanetions of crystal axis, temperature, and 
germanium concentration of the solvent. Plots of germanium dissolution 
versus time are generally characterized by two regions indicating an initial 
rate and a final rate. For the initial rates the [111] crystal direction is 
slowest and the [110] crystal direction is fastest. As the dissolution pro- 
ceeds, other crystal faces become exposed so that the final rates, beiyig less 
a. function of arty one crystal direction, are more nearly equal. When the 
dissolution rates are varied by temperature or by solvent constitution, it- 
is found that as the over-all dissolution rate increases, the rates in the 
three major crystal directions become more nearly equal. There^ are simi- 
larities between the dissolution of germanium by molten indiuni and the 
etching of germanium by an acid etch. 

Introduction 

MOLTEN indium pellets alloyed in the customary manner^ onto 
t^ermanium surfaces with low ed^^e dislocation density tend 
to spread excessively.^ Since bulk dissolution properties may 
be expected to play an important role in alloying, one of the approaches 
to coping with the problem of excessive spreading lay in modifying 
that part of the alloying procedure in which the germanium is dis- 
solved by the molten indium. Accordingly, experiments were performed 
to answer certain specific questions relating to such dissolution. These 
experiments include a study of dissolution as a function of crystal 
axis orientation, temperature, and amount of germanium already in 
the molten indium solvent. 



* Presented at The Electrochemical Society Meeting, Cleveland, Ohio, 
October 3, 1956. 

1 C. W. Mueller and N. H. Ditrick, ^'Uniform Planar Alloy Junctions 
for Germanium Transistors,'' RCA Revietv, Vol. XVII, p. 46, March, 1956. 

- JL V. Kettering and J. I. Pankove, "Crystal QuaHty as a Factor in 
Controlling the Alloying of Indium onto Germanium Surfaces," presented 
at the I.R.E.-A.I.E.E, conference on Semiconductor Research, Lafayette, 
Indiana, June 26, 1956. 
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Experimental Results 

Thin wafers of low-dislocation-density germanium about one half 
inch long, one eighth inch wide and about 40-60 mils thick were cut 
and carefully oriented so that the large-area faces were perpendicular 
to a major crystallogi^aphic axis. The wafers were surface ground and 
then etched in CP4 to remove about three mils of worked surface. (This 
was later increased to five mils; however, it was found that the results 
were independent of the amount removed by the CP4 beyond the orig- 
inal three mils.) The wafers were coated with a flux* to insure good 
uniform wetting and then partially dipped into the solvent (pure 
indium or a germanium-indium solution). Constant, nonviolent agita- 
tion was supplied to prevent the germanium already in solution from 
interfering with the dissolution process. After a certain time the 
wafers were removed, and the excess indium dissolved by HCl. No 
trace of recrystallized germanium could be detected. 

The dipped portion of the wafer was then compared to the undipped 
portion under a microscope. In the initial stages of the dissolution 
(5-10 seconds) a pattern is formed on the wafer surfaces consisting 
of plateaus of germanium which had been only slowly attacked by the 
indium, surrounded by regions of more rapidly attacked germanium. 
These plateaus vary in size from one micron to several mils. They are 
mostly circular with some few of the smaller ones triangular and 
hexagonal. With time the plateaus get smaller, finally disappear, and 
an over-all decrease in wafer thickness is observed. The surface at 
this point is not flat, and the decrease in thickness can be recorded 
only as a gross, averaged quantity. The data is plotted as the amount 
of germanium dissolved (in terms of the decrease in thickness) versus 
time. The measurements of decrease in thickness are accurate to about 
±10-2 centimeter. 

Figure 1 shows the rates of dissolution of germanium by pure 
indium at 347°C as functions of the primary crystal axes. The rates 
are given in Figure 1 and also in Table I. Each point on a curve 
represents a separate test wafer. These, and subsequent curves are 
characterized by two regions indicating an ^initial" and ''final" dis- 
solution rate. To compare initial and final rates, straight line segments 
have been used together with a "knee" to indicate the gradual transi- 
tion from the initial to the final region. 

A possible reason for the two regions in the curves may be thought 
at first to be due to an insufficient CP4 etch leaving some worked sur- 
face to be dissolved by the indium before the undisturbed bulk of the 



* Similar to the flux described by Mueller and Ditrick in Reference (1). 
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Fig-. 1 — Dissolution of germanium by molten indium, 7"=347°C. 

crystal was reached. However, when 5 X 10-3 centimeter, of ger- 
manium (the *'break point" in the curves of Figure 1) was etched in 
CP4 in addition to the amount ordinarily etched, essentially the same 
curves are produced. This ruled out the possibility mentioned above. 
It is now thought that the two regions generally observed in the dis- 
solution rates are due to the fact that although the test wafers have 
flat oriented surfaces initially, as the dissolution proceeds the sui^faces 
become more and more roughened. This means that ultimately the 
direction of the indium's attack is distributed among several crystal 
directions and not solely in any one direction. Thus, the initial regions 



Table I — Kates of Dissolution of Germanium by Molten Indium 





r r:z347°C 




Rate, 


in cm X 10-^^/min. 




Crystal 






Orientation 


Initial 


Final 


(110) 


10.0 


7.2 


(100) 


6.3 


5.0 


(111) 


3.1 


7.0 
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in the ciu'ves of Figure 1 are much more the property of crystal axis 
than the final region. It is found that, consistent with this, differences 
in the dissolution rates observed for the final regions are much less 
than differences in rates for the initial regions (see, for example. Table 
I). Note that in the initial rates direct experimental evidence is found 
for something which had heretofore been only strongly suspected — ■ 
namely that the formation of flat alloy junctions on a (111) crystal 
plane is due to the fact that indium dissolves germanium slowest in 
the [111] direction. 




TIME-SECONDS 

Fig. 2 — Dissolution of germanium by molten indium, [111] direction. 

The final rates given in Table I are little different from one another, 
the greatest difference being between the [111] and the [100] direc- 
tions, a factor of 1.4. However, the initial [110] rate is 3.2 times as 
fast as the initial [111] rate, and later on it is showTi that when the 
actual conditions of junction fabrication are more closely approximated, 
the ratio of the (110) alloying rate to the (111) alloying rate is even 
greater. 

The dissolution rate of the (111) oriented test wafers was then 
determined at 373°C, 417°C, and 498°C. This data is shown in Figure 
2. At 373°C and 417°C initial and final regions are again present, 
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with increased rates for both. At 498*'C only one rate appears although 
an initial rate may very well be present but not detectable by our 
methods because of the rapidity of the indium's attack. The dissolution 
rate increases sevenfold as we go from 347 °C to 498 °C. When the 
available data of Figure 2 is plotted as dissolution rate versus inverse 
absolute temperature, the results are as shown in Figure 3. Here, we 
see that the dissolution rate depends exponentially on the inverse abso- 




INITtAL RATES 



1.3 10= 14 



Fig. 3 — Temperature dependence of dissolution rate of germanium by 
molten indium, [111] direction. 



lute temperature. If we assume that the temperature dependence of 
the rate-determining step in the dissolution process is given by a Boltz- 
mann factor, then the activation energy associated with the disso- 
lution process is .53 electron volt. Interestingly enough, this is an 
activation energy associated w^ith the etching of germanium by No. 2 
acid etch found by P. Camp.^ Furthermore, on comparing the etch 
data of Camp's Table II with our owm data (Table I) it is found that 
the dependence of the dissolution rate on ci-ystal axis is the same 



3 P. J. Camp, "A Study of the Etching Kate of Single-Crystal Ger- 
manium," Jour. Electro, Chem. 5oc., Vol. 102, p, 586, October, 1955. 
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(qualitatively) for both No. 2 etch and molten indium. This behavior is 
true for both the initial and the final rates.* Thus, even though acid 
etching involves the formation of chemical compounds while dissolution 
by indium does not, identical activation energies and identical depend- 
ence of dissolution rate on crystal axis suggest that both acid etching 
and dissolution by molten indium may be very similar in mechanism. 

Increasing the temperature of the molten indium increased the 
rate of dissolution as expected. An experiment will now be described 
in which a slower dissolution rate was produced not by decreasing the 
temperature, but rather by adding germanium to the molten indium. 
In this way the conditions present during the fabrication of an alloy 
junction are more nearly approximated, since in the actual fabrication, 
which customarily involves very slow heating, the germanium is dis- 
solved by an indium-germanium solution. Accordingly, a mixture of 

Table II — Rates of Etching of Germanium by No. 2 Etch 



Rate, in cm. X 10-^/min. 



Crystal 
Orientation 


Initial 


Final 


(110) 


9.6 


6.0 


(100) 


8.4 


3.8 


(111) 


8.0 


4.3 



germanium and indium was prepared which produced, at the tempera- 
ture of the experiment, a solution about 75 per cent saturated with 
germanium. The dissolution rates were determined and the results 
are given in Figure 4. Essentially the same general characteristics as 
those found in Figure 1 appear. Two rates, an initial and a final one, 
are apparent in each of the curves. The final rates are less than the 
initial ones for the [110] and [100] directions, while the opposite holds 
for the [111] direction. Again we note that differences among the 
final rates are less than those among the initial rates. 

There are some important differences, however, between these 
curves and those in Figure 1. The general dissolution rates are about 
fifteen times slower than the rates using pure indium (note the abscissa 
scale change in Figure 4 from seconds to minutes). It is seen that 
with these much slower rates, the difference between the (110) rate 
and the (111) rate is given by a factor of five. Elaborating this par- 



* The identity for the initial rates of the two processes may not be 
significant. The initial rate with No. 2 acid etch may correspond to the 
removal of worked surface. In the case of dissolution by molten mdium, 
the worked surface had already been removed. 
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30 40 50 
TIME-MINUTES 



Fig. 4- 



-Dissolution of gei'tnanium by molten indium-gerinaniuni solution 
(75 per cent saturation). T = 347' C. 



licular point, the initial rates have been listed in Table III for those 
conditions where the over-all rates are mai'kedly different. The bottom 
line in this table gives the ratios of the (110) rate to the (111) rate. 
From this it is concluded that the slower the dissolution rate, whether 
it be due to temperature or to solvent composition, the greater the 
difference in initial dissolution rates, and consequently the greater 
the ratio of lateral dissolution speed to vertical dissolution speed on 
a (111) oriented surface. These are shown as the bulk dissolution 



Table III — Dissolution Kate of Germanium 



Crystal 
Orientation 


Rate, in 
In— 75% Ge 


cm X 10 

347 °C 
Indium 


~*V min. 

417°C 
Indium 


500 °C 
Indium 


(UO) 


.C)6 


10.0 


17 


--50 


(100) 


.50 




14 


- 50 


(111) 


.13 


3.1 


9 


- 50 


(110) rate 












5.0 


3.2 


1.9 


- 1 


(111) rate 
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processes B and C, Figure 5, in which the factors I'elating to the 
spreading of a pellet of molten indium during alloying are illustrated 
schematically. Process A, Figure 5, relates to wetting and surface 
tension forces. 



It has been found that surface undercutting, process B, Figure 5, 
is the major cause for excessive spreading of indium during alloying.* 
One way to minimize this effect would be to make processes B and C 
more nearly equal. The data in Table III indicates that this may be 
accomplished by increasing the over-all dissolution rate. This has, in 
fact, been found to be the case.'* When slow heating of the alloy units 
is used (as is customarily done) the germanium is dissolved by a 



Fig. 5 — Factors influencing spreading of a molten dot. Process (A) relates 
to wetting and surface tension forces. Processes (B) and (C) relate to 
bulk dissolution processes. 

solution already essentially saturated with germanium. This produces 
a very slow dissolution rate and a great disparity between the rates 
of B and C, Figure 5. When a very fast heating rate is employed, 
causing most of the alloying to take place at high temperature and 
with the dissolving solution far from germanium saturation, much 
faster dissolution rates and more nearly equal B and C rates are 
obtained. It is under these latter conditions that one finds no excessive 
spreading. It is to be emphasized that the fast heating rate per se 
does not produce the favorable alloying conditions, but rather the 
high temperature and nonsaturated condition of the solvent. We may 
remark in passing that even with the rapid heating, one can still 
produce the required flat junctions for which it was thought the slow 
heating was necessary. 



* G, W. Mueller, "Alloying Properties of Germanium Free of Edge 
Dislocations,** RCA Review, Vol. XVIII, p. 205, June, 1957. 
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ON THE NONLINEAR BEHAVIOR OF 
ELECTRON-BEAM DEVICES 

By 

F. Paschke 

RCA Laboratories, 
Princeton, N. .1. 

Abstract 

Summary — The Jioif linear space-char gc-wave equation is derived. By 
third-order successive approximation this complicated equation is split into 
three simultaneous linear differential equations, which are solved for the 
case of a velocity -7710 did a ted electron beam (klystron ). When excited at the 
fundamental frequency <d, the beam produces all harmonics with the fre- 
quency no). At fairly low levels the maximum current amplitudes are found 
to be proportional to the (1 • — -n/2)^^ power of the direct current. Equations 
for the gain and efficiency of a two-cavity klystron are given. The satura- 
tion efficiency is coynputed to he about 55 per cent. Thus the space-charge^ 
or ''longitudinal debunching/' does not reduce the efficiency of a klystron 
appreciably below the ell-known ballistic theoretic value of 58.2 per cent. 
Furthermore, it is shoivn that a computation of kinetic poiver floiv, using 
the linear solutions, is not quite justifiable. With the simultaneous linear 
differential equations given, it is possible to treat nonlinear phenomena, in 
traveling wave-tubes. 



T 



Introduction 

^HE theoretical evaluation of efficiency, harmonic frequency 
generation, mixing, etc. in electron-beam tubes requires an 
understanding of the nonlinear behavior of the device. This 
problem has been the subject of many papers, each dealing with one 
special type of tube, e.g., klystron and traveling-wave tube. Webstei-^*" 
investigated a velocity-modulated electron beam (klystron) under con- 
ditions of arbitrary signals, but took into account only the linear 
space-charge term. Hence, this theory is, for large signals, a purely 
ballistic one. Gittins and Sullivan,^ hoM^ever, showed experimentally 
that w^ith any appreciable current the neglect of space charge in a 



1 D. L. Webster, ''Cathode-Ray Bunching,'' Jom\ Appl. Phys., Vol. 10, 
p. 501, July 1939. 

2D. R. Hamilton, J. K. Knipp, and J. B. H. Kuper, Klystrons and 
Microwave Triodes, McGraw-Hill Book Company, Inc., New York, N. Y., 
1948. 

3 J. F. Gittins and A. B. J. Sullivan, "On Electron Bunching by Travel- 
ing Waves," S.E.R.L. Tech. Jour., Vol. 5, p. 127, December 1955. 
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velocity-modulated beam is never justifiable. Noi-dsieck/ in his travel- 
ing-wave-tube theory, followed a typical set of electrons and, again 
neglecting the space charge, computed their velocities and positions by 
numerical integration. Tien, Walker, and Wolontis^ and Rowe^ included 
space-charge effects by superimposing the space-charge force and the 
force of a divergence-less electric-circuit field. The space-charge field 
was calculated by replacing the delay line (helix) by a perfectly con- 
ducting drift tube which may lead, according to Reference (7), to 
appreciable errors. In addition, certain assumptions were made for 
the distribution of space charge along the axis. Tien, Walker, and 
Wolontis computed the space-charge force from a disk-charge model. 
Rowe presumed a sinusoidal space-charge variation along the axis. 
Neither presumption holds in actual practice. 

A straightforward way of calculating the nonlinear behavior of 
electron-beam devices is to take the electron stream not as a certain 
limited number of electrons but as a *'fiuid'* where the discrete charges 
are thought to be "smeared out" and to integrate the nonlinear space- 
charge-wave equation for the given boundary conditions. Here no 
superposition of a space-charge field and a circuit field is possible 
because of the nonlinearity of the field equation. In Maxwell's field 
theory there is only one definite total electrical field strength for which 
it is valid to say that 

div E = p/€q 

inside of the beam, and 

divjE'^iO 

outside of the beam. This "Eulerian" approach, however, is valid only 
until electron overtaking occurs which means that the velocity becomes 
a multi-valued function of space. As will be seen later, the theory 
always yields (for a given position) definite velocities. This means 
that in our fluid model overtaking does not occur. 

The nonlinear space-charge-Avave equation is very complicated and 



^ A. Nordsieck, "Theory of the Large Signal Behavior of Traveling- 
Wave Amplifiers," Proc, T.R.E,, Vol. 41, p. 630, May 1953. 

5 P. K. Tien, L. R. Walker, and V. M. Wolontis, '*A Large Signal 
Theory of Traveling-Wave Amplifiers," Proc, LR.E., Vol. 43, p. 200, March 
1955. 

8 J. E. Rowe, "A Large Signal Analysis of the Ti-aveling-Wave Ampli- 
fier," Trans, LR.E. PGED, Vol. 3, p. 39, January 1950. 

7 F. Paschke, **Die Wechselseitigkeit der Kopplung in Wanderfeld- 
rohren," Arch, Elektr. Uebertr., Vol. 11, p. 137, April, 1957. 
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extremely difficult to integrate. On the other hand we do not need a 
complete sohition, correct for all power levels. In order to predict 
saturation and the generation of the first two harmonics, a third-order 
theory is sufficient. The purpose of this paper is to derive the non- 
linear space-charge-wave equation and to solve it by third-order suc- 
cessive approximation for the special case of a drifting velocity- 
modulated beam (klystron). 

The Space-Charge-Wave Equation 
To simplify the analysis we make the following assumptions: 

1. One-dimensional confined electron fiow. 

2. The electric r-f field E in the direction of the stream is approxi- 
mately constant over the cross section of the beam. 

3. The electron velocity is small compared to the velocity of light. 

4. The effects of the potential depression caused by space-charge 
across the beam and the thermal velocity disti'ibution are 
negligible. 

The convection-current density 

io + t= (^0 + p) (Vo + v), (1) 
*o = Po^o 

appears in the continuity equation 

di dp 



The subscript 0 refers to the unpertui'bed d-c quantities at the 
input plane. The perturbations i, p, and v are caused by the r-f field 
and are functions of time and distance from the input plane. According 
to the assumptions listed above, the electrons are influenced only by 
the axial electric field E. Thus the force equation is 

+ (Vq + V) = yE, (3) 

dt dx 

The coupling between the beam and the structure is produced by 
the transverse electric field; hence, a displacement current per unit 
length, M, flows from the beam (Figure 1). This coupling current 
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enters into the continuity equation for the total current flowing through 
ail infinitesimal beam element; 



ZxM S 



(4) 



We introduce now a new^ variable, the "displacement'* 

y = €QE + m, (5) 
where m represents the transverse coupling field and is given by 

d-m M 



'dxdt 



(6) 




Fig. 1 — Current flow in an infinitesimal beam element. 

Equations (l)-(6) yield the following nonlinear differential equation 
for the displacement y: 



9-2/ 



Zy \ I Z^j 



-dy Y Zhj I -dy dy 

Po + — • + '^u 

Zx / 'dxU \ dx 



dt 



dxZt 



3-2/ \ / 

+ n • ( ^0 + 

dx- I \ 'dx 



Po + • 



Zy 



Zx 



dy 



dy \ I dy dy \ dhj y 

— — + '^^0 — — + — (y~m) + - — ) =0. 

dt \ dt dx I dx- Co \ dx 

a) 



The velocity v, i.e., the total velocity minus the initial d-c value, is 
related to the solution of Equation (7) by 
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dy dy 

+ ^^0 

dt "dx 

v= . (8) 

dy 

PQ-^ 

For calculations of efficiency only a third-order theory is required. 
Hence, terms of fourth and higher order in Equations (7) and (8) 
may be neglected without making appreciable errors if the discussion 
is restricted to powers below and up to saturation level. Under this 
assumption Equation (7) yields 

9-7/ 2 d^y 1 d^y 

+ + ^ + ^p2 (7j-m) 

dx- Vq 'dxM Vq"^ 

dy d-y dy Z-y \ 2 / dy d-y dy d~y 
+ 



o*- 0 



dx dt^ dt dxdt / \ dx dxdt dt dx"^ 



Sy dy 



€oVq^ dx 



d^y I dy \ 2 32y dy d^y ( dy 
+ — I — -2 . + 



V \ \ / dxdt dx dt dt^ \ dx 



Srj ( dy 



2 



+ (2/-m)=0. (9) 



Here /5p = o>^/Vo is the plasma-wave number. For very low power levels 
we can even neglect second-order terms. Under such a condition all 
terms in Equation (9) except the first four can be omitted; this leads 
to the well known space-charge-wave equation for small signals.^-^ We, 
however, try a third-order approximative solution 

y = Vi + ?/2 + 2/3 + * * • 

m = nil + 1712 + + • • * (10) 
v = -f + '^3 • 

Again omitting terms of fourth and higher order we get 



8 W. C. Hahn, *'Small Signal Theory of Velocity-Modulated Electron 
Beams,*' General Electric Revietv, Vol. 42, p. 258, June 1939. 

9 S. Ramo, *'Space-Charge and Field Waves in an Electron Beam," 
Fkys, Rev., Vol. 56, p. 276, August 1939. 
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+ (11) 

where 

'd^Vx 2 a^^i 1 Z^j^ 

/i= + + • + 

'dx^ Vq dxdt Vq^ dt^ 

/2 = + + + i^;- (2/2 " >»2) 

2 / dy^ dhj^ dy^ \ 2 / dy^ 3-*?/, 



(qVo \ 'dx M dxdt J \ dx 'dxdt 

'^Vi ^^2/1 \ 3?/ dyi 

+ (2/1 — ^%), 



3i dx^ I fQ^o" 
3% 2 92^3 1 327/3 

/s = V + — ■ ■ + (y-A ~ 'i^h) 

3a;2 z^zt ^2 3^2 

2 / 3'2/i 3^7/2 3?/o 3-?/i 37/2 ^^2/1 ^Z/i 3^7/^ 
+ . + _ 



IqVq \ 3a; 3^^ 3i- M 'dxdt dt dxdt 

2 I dy^ 327/2 37/2 3^7^^ ^y^ 32772 37/2 32?/^ 
+ + 



{q \ 3.T dxdt dx dxdt dt 3^2 2>t 3x'2 

1- / ^^2/1 / ^2/1 \- ^2/1 327/1 327J1 / 37/1 
- 2 -f 



Iq^ \ ?>x^ \ / 9^ 3a;3i 3^2 y ax 

37; / 3?/2 37^1 \ 37; / 37Ji 

(7/, — m-^) 4- 7/2 H • — Wi) 



To simplify the analysis considerably we set the terms of etiual order 
equal to zero. This procedure, of course, eliminates some of the pos- 
sible solutions of P"]quation (11). However, as will be shown later, the 
remaining solutions are sufficient to satisfy the boundary conditions 
in practical cases. The equations are 

/i==0, (12a) 

/2 = 0, (12b) 

/8 = 0. (12c) 
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Similarly we find from Equations (8) and (10) for the velocity terms 
Vi, and v-j^ the relations 



2i?/2 ^^/2 
■po^^2 = f-^'o ^'^'r 



(13a) 



(13b) 



^2/3 '^vs '^yi 

+ t'o + v. + ■ 



(13c) 



dx 



'dx 



dx 



Equations (12) represent a system of lineai' differential equations and 
the procedure of computing is now clear; we must first solve Equation 
(12a), put the solution into (12b), solve this Equation and put the 
solutions of (12a) and (12b) into Equation (12c) and solve it. The 
total solutions + ijo + Vs and 4- 1^2 + ^"^^^ satisfy the boundary 
conditions. This approach to the solution of a nonlinear differential 
equation is known as the successive-approximation method. 

Velocity Modulation by a Sinusoidal R-F Voltage 

Let us assume a cavity with very narrow gap at x = x^ modulating 
a beam of velocity ?'oo. The voltage across the gap is given by 



Vrrn—Fisin (o)i + ^). 



(14) 



From the conservation of energy we find, for the total velocity of the 
beam leaving the gap, 



dx 



1 sin (oit -(- (j)) . 



(15) 



dt /x^ V Voo^ 
Developing Equation (15) for the case 

2vV, 



into a series and omitting again terms of higher than third order, we 
get 
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sin (tD^ + ^) 
(16) 



+ 



cos 2 {U + 4>) + 



sin 3 (o)^ -F (/>) 



CURRKNTS AND VELOCITIES IN A MODULATED BEAM 

We assume a velocity-modulating cavity with a very narrow gap at 



(17) 



Equation (12) must be solved for m — 0, i.e., no coupling to the out- 
side and a velocity at a? = a;o given by Equation (16). After basic but 
rather involved algebraic manipulations, the solution for the displace- 
ment y is found to be 



— y = cos p^x sin {(Jit — PqX) 



V. 



00 



'^00^ 



(1 + cos 2 /3pX) sin 2 (oit — p,x) 



+ cos ppX cos 2 — ^^.t) 



-f ■ sin 2 PpX (1 — cos 2 ((at — /S^a:) ) 



16 \ V 



"jjO 



1 + 9 



■ cos 3 /3j,x sin (wi ~ /^^x) 



3 — 13 



' cos sin (ttii — 



+ 2 



sin 2 y^pO; sin (<oi — jS^x) 
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H (1 + cos 2 13 px) cos (isit — 13 ^x) 



3 (sin 8 f3pX + sin ^^o;) cos (od — ^c-'^) 



3 /vV, 



16 \W 



1 + 



(cos 3 f3pX + 3 cos /3pX) sin 3 (o)i — 



(sin 3 f3pX 4- sin l3pX) cos 3 (oi — I3^,x) 



+ (1 + cos 2 ppX) cos 3 (a — p^x) 



-\- - — [ j sin 2 PpX sin 3 (o>i — /^eX) 



cos sin 3 (od — f3^,x) 



The velocity, as derived from Equations (13), is given by 

1 fvViY vVi 

-f sin f3pX sin (o>t — ^^rr) 



(18) 



4 \ ?J 



1 fvViY o> r ^ 

— -sin2^, 



X sin 2 (o)i — ^^x) 



■ sin ^pO: cos 2 (coi — ^^aO 



+ (1 + cos 2 ^^a;) (1 — cos 2 (coi — ^,a;) ) 



] 



1 / r?yi \ 3 



16 \ Vnn^ / \ 



7 

1 / / '"pO 



sin 3 l3pX sin (o>i ~ I3^x) 
sin jSpO; sin (wi — /3^x) 
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4- ( 3 2 ( ■ ) j sin 2 p^x cos (wi — ji^x) — 2 

(1 + cos 2 ppX) sin (<ot — p^x) 



-\- I 5 2 I ■ ) cos 3 ppX cos {ttd — p^x) 



+ ( 3 + 2 1 1 cos p^x cos (o)i — 



sin jSj^a; sin 3 (o>^ — p^x) 
— ( 5 ■ 2 (cos 3 p^x + sin 2 ^^.r) cos 3 d^t — j^^.r) 



pO / ^'^pO \ ^ \ / ^'Vo 



~ 3 + 2 ■ cos Pj,x cos 3 Ut — p,x) + 2 



(l+cos2^pa:) sin 3 (o)i — J • (19) 

The boundary conditions at the excitation plane are met if </> in 
Equations (14)-(16) is set equal to 

<f> = . (20) 

2 0)^ 

The plasma-frequency oi^o in the amplitude terms of Equations (18) 
and (19) is the value before the beam enters the modulating cavity. 
Behind the cavity the plasma frequency is slightly changed and de- 
pends upon the amplitude because of the relation between d-c velocity 
and amplitude given by Equation (16). 

Beside the well known low-level waves^'^ with the propagation 
constants 



(21) 
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we have at large signals additional waves with the propagation con- 
stants 



<f) o) it: 2o)p at ± 3o)p 

for the fundamental frequency; 

2o) 2(1) ± (()„ 2(1) ± 2(1 



(22) 



7 = , , 

for the first harmonic; and 

3o) 3o) ± o>p 3r,> ± 2o,^ 3o) ± 3(-)^ 



(23) 



'^0 



(24) 



foi' the second harmonic. 

In most practical cases we may neglect high-order-space-charge 
terms in Equations (18) and (19) since 

«1, (25) 



and 

^v, Y 



(26) 



Under these assumptions we now compute the current density which 
is simply the negative time derivative of the displacement y: 

. — _ cos ppX COS (o)t — pgX) 



1 /vVA' 



(0 



+ . . (1 + COS 2 ppx) cos 2 (u>t — P,x) 

2 \ ^'oo- / \ ^Vo / 

1 / r,Fi \ 'V ^ \' 

+ ■ (cos 3 PpX + 3 cos p^,x) cos {id — p,x) 

32 \ too- / \ ^Vo / 



9 / 77V1 \ ^ 



^ 1 \ / 



32 \ W 



3 



(cos 3 PpX + 3 cos PpX) cos 3 {od — p,x) . 

(27) 
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By taking into account Expressions (25) and (26), Equation (19) 
yields the following simplified equation foi' the velocity: 

= sin p.pX sin {mi — p^x) 



V. 



00" 

1 / 7;yi 



4 Vt'oo^/ 0)^0 



7 / ^V, 



sin 2 Pj,x sin 2 {mi — p^x) 

(sin 3 ppX + sin p^x) sin (cot — p^x) 



2 



H I 1 I ) (sin 3 PpX + sin pj^x) sin 3 (a>^ — p^,x) , 

(28) 

Figure 2 is a plot of the absolute values of current density and velocity 
of the fundamental frequency versus distance for different power levels. 
The absolute maximum of current density (saturation) appears at 
x — ^ for a power level given by 

(29) 



^00 ^f>0 



The first set of curves corresponds to small signals, the third to satura- 
tion level. The fourth curve shows the amplitude distribution for 
higher than saturation level; the validity of this curve is doubtful 
because of the neglect of higher than third-order terms. Figure 3 
shows the amplitudes of the first two harmonics versus distance. Again 
the limitation to third- and lower-order term does not alloAV us to 
predict the saturation behavior of the harmonics. 

The question arises as to how large the power levels can be and 
still have Equations (27) and (28) represent good approximations of 
the exact solutions. To extrapolate to higher harmonics let us consider 
the **quasi-linear" case, which means neglect of saturation effects but 
not of the harmonic generation. The maximum current density (Equa- 
tion (27) ) appears at ct; = 0 and for quasi-linearity is given by 

ri\\ a> / ^ y 9 / <^ \® 
- — • cos toi -f ( I cos 2 tiit I ■ ' 1 cos 3 o>f . 



00 *"pO \ ^00 ^pO / ° \ ^00 ^'Vo 



(30) 
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Fig-. 2 — Ainplitudes of the fundamental current density (solid) and velocity 
(dashed) for different power levels versus distance. 

This equation seems to be the first part of the series 




Fig. 3 — Amplitudes of the harmonic current densities (solid) and velocities 
(dashed) versus distance. 
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Figure 4 shows this spectrum for different power levels. The series 
is convergent only if 



rjViii) 2 

<■ — = 0.736. 

2 „ e 



(32) 



00 



^pO 



Thus we may conclude that Equation (27) is a valid approximation 
for the total current density only for power levels well below the limit- 
ing value given by Equation (32) . This does not necessarily mean that 
the first terms of the power series for the f undamental current at a; — 0 
which can be obtained from Equation (27), 



1 



^00 ''Vo 



cos Uity 



(33) 

are applicable only for power levels given by Equation (32). This 
series seems to converge much more rapidly than Equation (31). At 
saturation level (Equation (29)), for example, the second term is 1/3 
of the first term. It should be noted that Equation (33) consists of 
the first two terms of the Bessel function of first order. This vague 
exti*apolation 



1 



COS o>£ ~ - 



192 



^pO 



(34) 

COS oit 



leads to a relation which is formally the same as the well-known equa- 
tion of the ballistic theory, i 



= 2Ji ( • ) cos ml. 



where z is the distance from the plane of excitation; 



(35) 



(36) 



This equivalence explains the good agreement between some experi- 
ments and ballistic theory.^ 
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To recapitulate, the total cim^ent density or velocity can be com- 
puted from Equations (18) and (19) but only for levels well below 
the limit given by Equation (32). We can, hoAvever, predict the be- 
havior of the current and velocity of the fundamental frequency up 
to saturation level. 

Comparison with the Ballistic Theory 

Webster^s ballistic theory^ predicts convection currents in a velocity- 
modulated beam which are given by 



> J„ • p,z COS n (od ^- (/>) . (37) 



i « / rh^V^ 

= 2 

in 



This result is valid for negligible space charge and our theory should 
agree in this case with Webster's. Shifting our coordinate system by 
the transformation of Equation (3G) we get from Equation (27) 

i 7)Vi oi 
_ . ^ . yii'i ^^2; cos (o)t - j3^x) 

^ [ . . j ^\]i2 p^^^ 2 (od — P^^X) 

1 fvV, 0, 

-\ 1 ■ ^ ) sin-^ jS^^: cos (oit — /3^x) 



8 \ r,r <'J;> 

7) Vi W \ 3 



1\~ O) 



^m^ f^^,z cos 3 (ud - - p^x) . 



(38) 



In Equation (38) the plasma frequencies and oj^q were considered 
equal because Equation (37) is based on the assumption of a pure 
sinusoidal velocity modulation with the d-c velocity unchanged. For 



— , (39) 
2 

the space-charge force certainly is very small since there is no bunch 
formed yet. In this case we get from Equation (38) 
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+ 



Pe^ 1 



COS 2 {id — jS^X) — ■ 



COS (id — P^X) 



(40) 



|8,,z cos 3 (wi — jSp.T) . 



The amplitudes of the currents of frequencies ??(.> are given by 



1: 



to 

H 



r,V, / 1 

■ p,z 1 

v,,^ \ 8 



■p,z 



(41a) 



(41b) 



(41c) 



From Equation (37) on the other hand, by developing the Bessel- 
functions into a power series for small arguments, 



2l>n 



nl 



/ nr)Vj 



n + 1 



(42) 



As expected, the Equations (41) and (42) are in agreement. Ac- 
cording to measurennents of Gittins and Sullivan^ the neglect of space 
charge with any appreciable beam current is not justifiable except for 
positions close to the excitation plane where Equation (89) holds. 

Two-Cavity Klystron 

We assume a narrow-gap cavity at a*— — 7r/(2^p) which produces a 
velocity modulation given by Equation (10). At the output cavity, 
which is assumed to be at a; = 0, the current density of the fundamental 
frequency is given by Equation (33). We now introduce the beam 
voltage 



Vim'' 



2v 



(43) 



Saturation of the current density is reached at 
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Vi sat <^ 

where 



1/ 



(44) 



2 \ 

^2z-o( J . (45) 

In order to prevent reflection of electrons, the Q of the second cavity 
must be adjusted so that the r-f gap-voltage does not exceed the d-c 
beam voltage. For maximum output the gap voltage equals the d-c 
voltage, and the saturation power delivered from the beam to the 
cavity is given by 

1 

/^2sat-— \h.at\ SVo. (46) 

2 

Here S is the area of the cross section of the beam. The electronic 
saturation efficiency is given by 

^2 sat 1 I ^1 sat I 
Ve\ = ~ — ■ ] • 

From Equation (45), 



2 \3/2 

V.i = [ " =^54.4%, (47) 



If we use the vaguely extrapolated Equation (34), the calculated satu- 
ration efticiency agrees with the result of the ballistic theory, i.e., 



7?,i = 58.2%. (48) 

In any case we can conclude that the space charge or the so called 
'iongitudinal debunching" has no appreciable influence on the electronic 
efficiency. This is in agreement with experimental results.^ For other 
data, however, such as velocity- and current-modulation or gain, the 
ballistic theory only holds near the plane of excitation as previously 
discussed. 

If the equivalent gap resistances Ri and R2 of the input and output 
cavities are known, the gain of the klystron can be calculated. The 
input power, P^, is given by 
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= , (49) 

and the output power by 

A'.. (50) 

2 

Equations (83), (48), (49), and (50) yield an expression foi* the 
gain, 




(51) 



where /n is the d-c beam current, 

h = (52) 
Thus the small-signal gain is given by 




and the saturation gain (from Equation (44)) by 




(54) 



From Equations (53) and (54) one can see that the ratio of saturation 
gain to low-signal gain is 4/9, which corresponds to a difference of 
3.52 decibels. 

Because of the limitation of Equation (32), the harmonic output 
can be predicted only for levels far below saturation. Figure 4 shows 
the spectrum for three different power levels which is given by Equa- 
tion (31). While in the ballistic theory all the harmonic-current ampli- 
tudes depend linearly upon the direct current (Equation (37)), in 
the theory presented here a relation 
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n 

1 

is predicted. The fundamental current is proportional to the square 
root of the direct current, as is well known. The first harmonic output 
is predicted to be independent of the direct current whereas the higher 
harmonic amplitudes should decrease with increasing direct current. 
This is understandable because the smaller the d-c beam current for 
a given r-f power level the more nonlinear the system becomes. The 
electron stream acts like a linear system if the r-f current is very 
much smaller than the direct current. 

Kinetic Power Flow 

The third-order solution makes it possible to clarify the question 
of kinetic-power flow in electron streams. This concept was the subject 
of a number of papers. ^""^-^ According to References (14) and (12), 
it is a consequence of the conservation of energy that 

dxvTj^^lE, (56) 

where P^r '^^ the time average of the kinetic-power flow P^^ which, in 
the most general form,'^ is given by 

1 

Pa' = C?o + 0 (^;o + ^')-. (57) 

Using Equations (18), (19), (43), and (57), 
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(58) 



In this case the d-c term 



is of no interest and is therefore omitted. Equation (56) can be 
checked by computing iE from Equations (4) and (5) (m = 0, il/ — 0) ; 




Fi^. 4 — Frequency spectrum of the convection-cui'rent density 
in a quasi-hnear electron beam. 



iE — — y — 
^0 



1 3(7/^) 



■V' 



(59) 



^0 



This equation yields a nonzero value for iE only if there is a d-c 
term of the displacement, hence, a d-c term of the electric field. 
Indeed, we find from Equations (18) and (59) 



iE=: —I ■ ^,sin2^^^. 



(60) 
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From Equation (58), 

div Pj, = — - ^ — ( j sin 2/3,x, (61) 



dx 8 \ F, 



which agrees with Equations (56) and (60). The time-average value 
of iE in Equation (60) is caused by space charge, and originates from 
a second-order term of the displacement y. Thus the calculation of 
kinetic power flow by use of only the linear solution is 7iot justifiable, 
even at arbitrarily low levels. For ion neutralization a d-c term of 
the electric field cannot exist. In such a case we have to omit the d-c 
term in Equation (18) and the kinetic power flow is zero. These re- 
sults are in agreement with those of Walker.^^ 

The physical explanation for the maximum of kinetic-power flow 



FAST ELECTRON 
DECELERATED BY 
SPACE CHARGE 

\ 




SLOW ELECTRON 
ACCELERATED BY 
SPACE CHARGE 



PLANE OF EXCITATION 

Fig. 5 — Electron-time table. 



at x = 0 is the increased average velocity of the beam at that point 
due to space-charge effects. In the '^electron-time table" (Figure 5) 
the ballistic paths of three representative electrons are plotted, one 
passing the excitation gap at maximum r-f field (fast electron), one 
at zero, and one at minimum r-f field (slow electron). Because of the 
higher charge density, the space-charge forces are obviously much 
more efficient for the slow electrons. Hence, the slow electrons are 
accelerated by the space charge more than the fast electrons are 
decelerated. The net effect in the equilibrium state of the plasma 
oscillation at a: = 0 must be an increase in average velocity. Indeed, 
Equation (19) yields the average velocity 

'^00+'^ 1 / y^V, Y 

■ =^ 1 + — — cos 2l3,,r, (62) 

"^00 4 \ VQ^y^ j 
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which has a maximum at a: = 0. The value at excitation plane 
(/3pX — ~ it/2) agrees, of course, with the boundary condition given 
by Equation (16). 

NOTE ON THE NONLINEAR BEHAVIOR OF TRAVELING-WAVE TUBES 

The system of linear differential Equations (12) allows us to 
calculate also the nonlinear behavior of traveling-wave tubes. 

If we make the simplifying assumption that at the harmonic fre- 
quencies there is no coupling between beam and delay line, it is obvious 
from Equations (12) that the nonlinear behavior of traveling-wave 
tubes can be described by the small-signal parameters, since the small- 
signal solution determines also the higher-order solutions and V:^. 

This seems to be in contradiction to the theory of Tien, Walker, 
and Wolontis^ who found that, at least for high levels, the plasma 
frequency reduction factor R must be added as a new parameter. It 
has, however, been showm that this parameter also appears at low 
levels because it determines the coupling of the space-charge wave to 
the delay line.^ 

Acknowledgment 

The author is greatly indebted to S. Bloom for reviewing the 
manuscript of this paper. 



www.americanradiohistorv.com 



A SURVEY OF METHODS USED TO DETERMINE 
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Sumntary — This paper discusses the nature of contact potential and 
the effects of changes in grid-cathode contact potential on "bias skiff' in 
receiving tubes. It describes the method of determining true grid-cathode 
contact potential by extrapolation of portions of the grid'Characteristic 
curve, ^ as well as two direct methods of measurement tvhich are generally 
used m preference to the extrapolation rnethod because of their relative 
simplicity. Although these direct methods do not indicate true contact 
potential, they can, in many cases, provide information of considerable 
value in bias-shift studies. The relationships among the three methods 
and the application for which each is best suited are discussed 

Introduction 

A PROBLEM which frequently confronts manufacturers and 
large-scale users of receiving tubes is a variation in the 
grid-bias requirements of an amplifier type. This variation 
(called bias shift) may occur from lot to lot, among tubes in a particu- 
lar lot, or during the life of an individual tube. Although normally of 
small magnitude, this variation can cause a substantial change in a 
major characteristic, particularly when it occurs in a high-mu triode 
or a high-transconductance type. It is particularly objectionable to the 
equipment manufacturer because it may necessitate frequent modifica- 
tion of circuit design to change the value of an applied bias, or the 
use of a compromise value of bias which limits the performance ob- 
tainable from the equipment. 

It is more or less customary to attribute bias shifts to changes in 
contact potential (the difference in potential developed between dis- 
similar metals which are in direct mechanical or electrical contact). 
This difference of potential is created by, and equal to, the difference 
in the work functions of the metals, and its polarity is such that the 
metal having the higher work function is negative with respect to 
the other. In the case of a tube, the electrode having the lowest w^ork 
function is usually the cathode. Consequently, the contact potential 
between the cathode and grid of a triode amplifier tube is equivalent 
to an applied negative grid bias, and changes in the work function of 
either the cathode or the grid wnll cause corresponding changes in the 
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effective grid bias. Such changes in work function may result from 
manufacturing variations in the composition of the grid and cathode 
materials, or from contamination of the grid by the cathode during 
the tube manufacturing process or the life of the tube. 

The evaluation of true contact potentials in tubes is complicated 
by the fact that the resulting bias voltages can be affected by the pres- 
ence of negative grid currents due to gas, grid emission, and leakage. 
Because there is no method of direct measurement which will eliminate 
the effects of all of these contributory factors, true contact potential 
cannot be measured in terms of the equivalent bias voltage. Neverthe- 
less, a number of methods for the determination of ^'contact potentiar' 
(so-called) have been developed and are used throughout the industry. 
Some of these methods ai^e based on direct measurement of grid-bias 
voltages or grid currents and, therefore, include the effects not only 
of true contact potential but also of potentials due to other factors. 
This paper describes the various methods now in use and their advan- 
tages, shortcomings, and proper applications. In the interest of 
standardization, it also suggests a system of nomenclature and sym- 
bolization which indicates the character of the results obtained by each 
of the methods described and the conditions under which the results 
are obtained. 

Work Function 

Before proceeding with a discussion of contact potential, it is 
desirable to review the concept of work function. An electron a very 
short distance from the surface of a conductor in a vacuum is sub- 
jected to a combination of image and quantum forces which tend to 
prevent its escape from the parent metal. ^"•-^'•^^'^ Minor quantum effects 
(such as ''tunnel" processes) Avhich permit passage of a few electrons 
through the work-function barrier will be neglected here. The elec- 
tron may pass through the surface of the metal, but it will return 
unless it possesses enough energy to overcome the restraining forces. 
This behavior is comparable to that of a ball which is propelled up an 
incline, and which will roll back if it has not been pushed hard enough 
to make it go over the top. When a metal is heated, the internal elec^ 
trons partake of the thermal energy in varying amounts, some getting 

1 G. Herrmann and S. Wagener, The Oxide-Coated Cathode, Vol. II, 
Chapman & Hall Ltd., London, England, 1951: (a) p. 13; (b) p. 1G8; 
(c) p. 192; (d) pp. 32, 188; (e) p. 190; (f) p. 34. 

2E. L. Chaffee, Theory of Thermionic Vacuum Tubes, McGraw-Hill, 
N. Y., 1933: (a) p. 56; (b) p. 80. 

3W. B. Dow, Fundamentals of Engineering Electronics, John Wiley 
Inc., N. Y., 1952: (a) p. 218; (b) p. 203. 
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enouRh to permit them to escape into the vacuum. The work function 
of a metal (the equivalent of the height of the incline) is the total 
energy which an electron at the Fermi level must possess to escape 
completely and is usually expressed in electron volts. A plot of the 
distance out from the surface to which electrons having given amounts 
of initial energy will go is called an "energy profile'' diagram. Such 
a diagi'am is shown in Figure 1. 

The energy profile representing the work function of an emitter 
is altered by the presence of an applied electric field between the 
emitter and an adjacent electrode. A difference in work function 
alone between an emitter and an electrode with which it is in electrical 
contact is equivalent to an applied field. Under such conditions the 




— EMITTER 



— DISTANCE — 

1 — Knorgy-profile diagram of a thermionic emitter in vacuum, in the 
absence of an applied field. Only electrons having energy levels higher than 
that indicated by the dashed line can escape; all others return to the emitter. 



energy profile may not have the form shown in Figure 1, but climbs 
or dips in accordance with the sign of the applied field. These effects 
are illustrated in Figure 2. In Figure 2a, curve I shows the retarding 
field produced when the emitter is in electrical contact with a collector 
electrode having a larger work function. Curves II and III indicate 
the effect of progressively increasing anode currents which super- 
impose a space-charge field on the retarding field. In curve II the 
additional effect of the space-charge field is not sufficient to interfere 
with the flow of electrons to the collector. In curve III, however, the 
current is assumed high enough to produce a space-charge field suffi- 
cient to limit the electron flow. Figures 2b and 2c show the effects 
on the energy profile of applied potentials having the polarities indi- 
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cated. Curves II and III in both figures show the effect of the presence 
of an increasing number of electrons in the space between the elec- 
trodes, curve III, in each case, showing space-charge limitation. Note 
that if the collector is sufficiently positive with respect to the emitter 
so that the current is limited only by the emitter work function, the 
resulting field is termed a ''saturating'* or ''temperature-limiting" field 
(Figure 2c, curves I and II). 




Fig. 2 — Typical energy-profile diagrams of a thermionic cathode in the 
presence of a collector electrode having a higher work function, showing 
the effects of various applied voltages. 

What has been said thus far applies to pure metal emitters but not 
necessarily to oxide-coated cathodes which are the types most commonly 
used in receiving tubes. The mechanism of emission from the oxide 
cathode is more complex than that from a simple metal and is not yet 
fully understood.!^ The work function of an oxide is a function of 
current^^ and is believed to depend on temperature as well.^^-^'' For 
most measurements on commercial electron tubes, however, the cur- 

4 L. S. Nergaard, "Studies of the Oxide Cathode/' RCA Revietv, Vol. 
13, December, 1952: (a) pp. 512, 520; (b) p. 496. 
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rent and temperature effects can be considered second-order effects. 
In this discussion of contact potential it is assumed that the simplified 
picture of work function applies. It is also assumed that the contact 
potential of interest is that between the cathode and the adjacent 
electrode. 

Methods of Measurement 

The methods used to determine contact potential fall in three 
general categories. In one, a potential is determined graphically by 
extrapolation of two well defined regions of the grid-current-grid- 
voltage characteristic curve. This procedure yields what will be con- 
sidered true contact potential. In the second method, a potential is 
determined from the applied voltage required to produce a chosen 
value of collector current. In the third, a potential is determined 
from the reading on a high-impedance voltmeter connected between 
grid and cathode. 

The oxide cathode, being a semiconductor, exhibits a thermovoltaic 
effect which modifies the contact potential and must be corrected for 
in exacting work. The magnitude of this correction is of the order of 
0.003 volt per degree of temperature differential across the coating.^^ 
There is also a small effective voltage due to the initial velocity of the 
electrons.^'' 

The three methods and their significance are described and an 
appropriate title and symbol suggested for each. 
Method 1 — Intersection Method 

In this method the cathode is heated to a temperature in the 
neighborhood of 600°-700^K by operation of the heater at 24 per cent 
of its rated voltage (for example, 1.5 volts for a 6.3-volt heater). To 
minimize the effect of the plate (assuming a triode), the plate is tied 
to the cathode. The logarithm of the grid current is then plotted as 
a function of the applied grid voltage. The resulting grid-characteristic 
curve usually contains three clearly defined regions representing the 
effects of the retarding, space-charge, and saturating fields. 

The technique by which such a curve may be used to determine 
true contact potential is illustrated in Figure 3, Assume that the 
cathode has a work function of 1 volt. Figure 3a gives the energy- 
versus-distance configuration corresponding to an applied grid voltage 
of —0.25 volt, resulting in an assumed grid current of 0.001 micro- 
ampere. Even though the grid current is quite small the electrons 
form a space-charge field, with the result that the potential distribu- 
tion between the grid and cathode is not linear. The space-charge field, 
however, does not form a barrier higher than that of the grid work 
function, so that the latter plus the applied field are the only factors 
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limiting the current. Only electrons having thermal energies greater 
than 1.25 electron volts can reach the grid to make up the 0.001 
microampere of current. 

Figure 3b shows the profile for zero applied voltage, i.e., the Fermi 
levels of the cathode and grid have the same energy value. Under 
these conditions the only retarding field affecting the electrons is that 
produced by the difference in work function, i.e., by the true contact 
potential. The grid, therefore, accepts some lower-energy electrons, 
with the result that its current increases. When the grid is made 




Fig. 3 — Diagram showing the principal regions of a grid-characteristic 
curve, and determination of true contact potential by the intersection method. 

positive with respect to the cathode (Figure 3c), the profile approaches 
the horizontal, and at a certain point the space-charge field begins to 
become important in limiting the current. Up to this point the grid 
characteristic plot is a straight line, because in the retarding-potential 
region the collector current from an oxide cathode, like that from a 
pure metal emitter, is generally a logarithmic function of applied 
voltage as predicted by theory.^^* The interposition of the space-charge 
barrier at this point, however, tends to reduce the current, with the 
result that the curve starts to round off (Figure 3d). 
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Fi4,nire 3e shows that if the grid is made sufficiently positive, both 
the grid work function and the space-charge barrier no longer have a 
limiting effect on the current and all of the electrons emitted are 
collected. However, the current in this saturation or temperature- 
limited region of the grid characteristic curve does not remain con- 
stant, but increases slightly with applied voltage due to a slight change 
of the cathode work function by the applied field. 

If no space-charge barrier were present, the retarding-field region 
of the curve would change abruptly into the saturation region at the 
applied grid voltage which sets the upper rims of the cathode and grid 
work functions at the same level (Figure 3d). This applied voltage is 
numerically equal to the difference in work functions and, therefore, 
is the true contact potential which we are trying to determine. Be- 
cause there is always a rounding of the curve due to space charge, 
this voltage, Ecp, is given by the intersection of the extrapolations 
of the retarding and saturation x^egions of the curve. In circuit anal- 
ysis Ecp must be preceded by a minus sign because it acts in the same 
manner and direction as a negative grid bias. 

Although the intersection method is the best one available for 
determining true contact potential in commercial receiving tubes, it 
has several shortcomings. (1) To avoid damage by saturation currents 
the curve must be taken at low cathode temperature. Because of the 
dependence of the cathode work function on temperature, it is to be 
expected that the contact potential would be different at operating 
temperature. (2) Resistive layers on the cathode base metal and on 
the grid may aflTect the distribution of the applied voltage, thereby 
distorting the grid-current curve.^ Consequently, even the retarding- 
field region of the curve may be rounded, making extrapolation difficult 
and leading to erroneous results. The curve can also be rounded by 
nonuniformity of the cathode work function, as evidenced by variations 
in emission density over the surface of the oxide coating.^« Some 
authorities believe that the shape of the curve may also be influenced 
by the porosity and roughness of the oxide surface. Fortunately, these 
adverse factors are not significant in most tubes and the values ob- 
tained when precautions are taken to assure clean electrode surfaces 
show good agreement with the expected contact potentials. (3) The 
intersection method yields true contact potential only for plane-parallel 
electrodes since the slope of the retarding portion of the curve varies 
with the electrode geometry.^^' ^'^ Fortunately, most tubes approach 



5 G. C. Dalman, Effects of Cathode and Anode Resistance on the 
Retardmg Potential Characteristics of Diodes/' Jour. Ann! Phvs Vol 2^ 
p. 1263, October, 1954. * *' 
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the plane-parallel configuration closely enough so that, to a first approx- 
imation, this consideration may be ignored. 

Method 2 — Potential Producing an Arbitrary Reference Current 

An applied potential which produces a certain chosen reference 
current in the retarding-field region of the grid-characteristic curve 
is called a retarding potential, and is symbolized E,,c, the subscript C 
indicating a relatively cool cathode (one operated at 24 per cent of 
the rated heater voltage). The reference current chosen should be low 
enough to assure that the corresponding applied potential is truly 
within the retarding-field region. A suitable value is 0.1 microampere. 
From Figure 3 it can be seen that the retarding potential which gives 
-fO.l microampere is determined only by the magnitude of the grid 
work function, the cathode work function becoming significant only in 
the saturation region. Differences in values of E;^,^, therefore, reflect 
differences in grid work functions, and are equal to differences in true 
contact potential only if the cathode work function can be assumed to 
be constant. 

If the +0.1 microampere reference current is outside the true 
retai'ding-field region because of poor cathode activation or large 
voltage drops across electrode resistances, differences in Eifc may not 
truly indicate differences in grid work functions. Consequently, the 
smaller the reference current that can reliably be measured (assuming 
no leakage currents are present), the more accurately will changes in 
retarding potential indicate changes or variations in grid work func- 
tion. It is especially important that measurements of retarding poten- 
tial made for comparative purposes be made at the same cathode 
temperature. Because increases in cathode temperature raise all por- 
tions of the grid-current-grid-voltage curve along the current axis, 
they have little effect on the contact potential but substantial effects 
on the voltage necessary to produce a given grid current. An increase 
in cathode temperature, therefore, results in a reduction of the value 
of Ei^c^ i-e., ^L'c becomes less positive, or more negative. 

Because retarding potentials lie well below the saturation region, 
they may be measured at normal cathode temperatures. Under these 
conditions it is necessary to apply a negative voltage to the grid to 
retard the majority of the available electrons. Figure 4 shows the 
location of the retarding potential point on a conventional grid-current- 
grid-voltage curve, the subscript H indicating a '*hot cathode" (normal 
temperature). £7,^^ will not only vary (like E^^r) with changes in grid 
work function but will also shift slightly with changes in grid emission 
and leakage currents, which increase with cathode temperature. It is. 



www.amerinanradiohistorv.nom 



CONTACT POTENTIALS IN RECEIVING TUBES 251 



therefore, not as accurate as Ej^.q as an indicator of contact potential, 
but it is, nevertheless, useful as a means of showing gross shifts in 
grid work function. 

A grid-current-grid-voltage curve taken with voltage applied to 
the plate differs from one taken with the plate floating or grounded 
because of negative currents resulting from gas ions and increased 
leakage currents' Retarding potentials measured under these condi- 
tions provide a general indication of bias shifts due to all causes, and 
are, therefore, only partially dependent upon contact potential. To 
avoid confusion with the retarding potentials Ej^fj and Ej^ji, which are 



Fig. 4 — Curves showing differences in limiting and developed potentials 
when tube is operated with plate floating or grounded (subscript E) , and 
with plate at an applied positive potential with respect to the cathode 

(subscript iV) . 

measured with the plate floating or grounded, a retarding potential 
measured with the plate positive is termed ''limiting potential'' and 
is symbolized as £'l(am> the subscript .Y representing the reference 
current in microamperes. Contact potential is sometimes defined as 
the potential required to limit the grid current to zero microamperes, 
i.e., ^//,(o)- Another value frequently used to indicate or determine bias 
shifts is E^,,o.n which has practically the same significance as £'/,(o)- If 
negative grid currents are used for reference, the limiting potential 
may be difficult to determine because of the instability of the negative 
component due to gas. It is undesirable, therefore, to reference the 
limiting potential to currents in the negative region. 




0.1 
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Method 3 — Voltmeter Method 

If a d-c voltmeter having a resistance of at least 10 megohms (such 
as a vacuum-tube voltmeter) is connected between the grid and cathode, 
and no other connection is made to the grid, a voltage drop will be 
developed across the meter by the grid current of the tube. Graphically, 
this voltage drop (called ^'developed potential") is determined by the 
intersection of the load line with the grid-current curve and, therefore, 
will have one value when measured with the plate at a positive potential, 
and another when measured with the plate floating or grounded. These 
two developed potentials are symbolized, respectively, as £7/>v indi- 
cating normal plate voltage) and Ej^y (F indicating a floating or 
grounded plate). It can be seen in Figure 4 that Ej^j,. and Ej^y differ 
very little from retarding and limiting potentials measured under the 
same conditions and, consequently, provide similar information. The 
developed potentials, however, are more easily measured than the 
others. 

Relationship Between Methods 

Because all three methods described above, with variations, are or 
have been used to determine contact potential, it is desirable to estab- 
lish as precisely as possible the relationships among them. The writer 
has participated in conferences of tube-industry scientific and engi- 
neering personnel at which at least three of the various potentials 
defined above have been repeatedly but inaccurately referred to as 
"contact potential". A typical example of this misuse occurred in a 
discussion of a manufacturing variation which resulted in a decrease 
in the grid work function, and a corresponding increase in developed 
bias. This increase was erroneously described as an increase in contact 
potential. The error of this description can be seen by reference to 
Figure 3, which shows that a decrease in grid work function moves the 
entire grid-characteristic curve to the left, and, therefore, reduces the 
contact potential. 

A decrease in cathode work function raises the saturation region 
of the curve, but does not shift the retarding-field line. An increase 
in cathode temperature enlarges the entire curve, particularly the 
rounded portion. Theoretically, cathode temperature also determines 
the slopes of the retarding and saturation regions, but in practice these 
slopes may also depend upon other factors.^' 

The curves of Figure 4 show true cathode-to-grid retarding-field 
currents combined with negative grid currents. A decrease in grid 
work function shifts the cuiwes to higher negative voltages and in- 
creases the absolute value of the developed, limiting, and retarding 
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potentials. An increase in ion currents due to an increase in tube gas 
content, however, makes Er.^x) ^^'^<^ ^ny Jess negative but does not 
affect the other measurements. Shifts in cathode work function have 
little effect on the values shown in Figure 4 since, at normal cathode 
operating temperatures, the primary grid current consists of electrons 
having energies higher than the energy of the cathode work-function 
barrier. 

Table T gives the results of measurements made on an experimental 
twin triode similar to type 6CG7 by each of the methods described 
above. The grids of the two triode units in the experimental tube 
were made of different materials, but were equally contaminated from 
the cathode coating. The differences, not only in magnitude but also 



Table 


I — Measurements on an 


Experimental Twin 


Triode 


Potential 


Triode Unit I (S 
(Magno-Nickel Grid) 
(volts) 


Triode Unit 2 
ilver-Plated Magno- 
Nickel Grid) 
(volts) 


Difference 
(Unit 2-Unit 
(volts) 


Eri' 


1.27 


0.50 


—0.77 


Erc 


0.92 


0.16 


— 0.7G 


JIj itn 


—0.06 


—1.35 


—0.39 


Euo.i) 


—0.91 


—1.13 


—0.22 


E Li(i\ 


—1.00 


—1.25 


—0.25 


Enr 


—0.94 


—1.31 


—0.37 


E},.v 


—0.85 


—1.13 


—0.28 



in the direction of change obtained with the various methods, emphasize 
the necessity for care in the choice of method used to determine bias 
shifts. 

Conclusions 

Each of the methods described above has a specific field of applica- 
tion in bias-shift studies. The intersection method is the only one 
which will indicate shifts in true contact potential, i.e., those resulting 
from changes in the work functions of the cathode, the grid, or both. 
The retarding-potential method is best suited for measurements of 
shifts due to changes in grid characteristics — Ei>(. showing the effects 
of grid contamination, and E/^j^ the combined effects of grid contami- 
nation, leakage, and grid emission. Limiting potential and developed 
potential show shifts due to all causes. To avoid confusion, the method 
used for obtaining bias-shift data should always be clearly identified. 
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Suntmary—Oxidafiov of the germanium staface is shown to be a major 
factor mfinenchig surface recombination and thereby the current transfer 
ratio, oich, of germanium ]}-n-p transistors. When the transistor is etchi^d 
electrolytically in. KOH, the surface becomes covered by a Jiydrated ge^'- 
maimnn monoxide-dioxide layer giving loiv surface recombination and high 
oLcb. DuHng the subsequemt life of the transistor this layer grows through 
further oxidation , the higher the temperature the more rapid the groivth, 
causing a reduction in ccrh. As the oxidation is rapid initially and then slows 
down ivith time, it has proven advantageous in manufacturing to subject 
netv transistors to heat treatment. This creates a certain thickness of oxide 
tvhich during the subsequent life of the transistor doei4 not change rapidly 
even at elevated temperatures. 



RIFT of the electrical characteristics of transistors, especially 



of the current-transfer ratio, a^f), has lowg been a major prob- 



^^^-^ lem in transistor applications. It has been known for some 
time that this drift is associated with changes in surface conditions.^ 
It is shown here that oxidation of the germanium surface is a major 
factor in the slow decay of a^^^ in p-n-p transistors. The rate of decay 
increases with temperature, and a^j, may change to such an extent that 
transistor circuits designed around the initial value of a^h fail to 
operate properly after a period of time. Our knowledge of the oxidation 
of germanium is still leather incomplete, but wuth the help of recent 
publications, an outline of the processes responsible for this behavior 
can now be sketched. 

Evidence is presented to support the picture that n-type germanium 
electrolytically etched in KOH has a low surface recombination velocity", 
a high and is covered by a thin GeO layer. The surface layer 
oxidizes further at temperatures between 25^ and 125°C according to 
the so-called **cubic" law, as is the case with copper. This oxidation is 
in the form of a GeOo layer M^hose growth is responsible for an increase 
in the surface recombination velocity and thereby the decay in the 
current-transfer ratio during life. 

^ R. H. Kingston, "Review of Germanium Surface Phenomena/' Jour, 
AppL Phys., Vol. 27, p. 101, February, 1956. 
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Earlier Work on Atomically Clean Surfaces 

On a freshly cleaved germanium surface, Green and Kafalas- have 
shown that at temperatures between —68° and +25°C, chemisorption 
in less than five seconds covers the surface with about three oxygen 
atoms per surface germanium atom, thus building up a very thin initial 
oxide layer. From then on the further chemisorption follows the usual 
logarithmic law,^ 

ATr^CilogGO^ + (1) 

where N is the number of oxygen atoms absorbed per cm^, 

t is time in hours, 
Ci = 1.4 X 1014, 
C2 = 8.4 X 1014. 

A logarithmic law^ was found for both p-type and n-type germanium, 
and could be followed for 22 days until the rate was too small to be 
observed. 

It has been shown by Law and Meigs^ that at 500 °C oxidation 
follows the parabolic law: 

iV ^ 60C3^i/2, (2) 

where C.-i 9 X lO^^ x v^^', V is the oxygen pressure, and t is the time 
in hours. This is typical of a process where diffusion through the 
formed oxide layer is the rate-determining factor. At still higher 
temperatures, Law and Meigs^ and Bernstein and Cubicciotti^ report 
a linear law in which the rate-determining factor is the evaporation 
of germanium monoxide from the surface: 

^ C^t, (3) 
where C4 is a pressure-dependent factor. 



2M. Green and J. A. Kafalas, ^'Oxidation of Glean Surfaces of Ger- 
manium Below 2'5°C," Phys, Rev., Vol. 98, p. 1566, 1955. 

3 M. Green, J. A. Kafalas, and P. H. Robinson, "Interaction of Oxygen 
with Clean Germanium Surfaces: Part I Experiment," Semiconductor Sur- 
face Physics, University of Pa., 1957, p. 349. 

T. Law and P. S. Meigs, "The High Temperature Oxidation of 
Germanium," Semiconductor Surface Physics, University of Pa., 1957, p. 383. 

5R, B. Bernstein and D. Cubicciotti, "The Kinetics of the Reaction of 
Germanium and Oxygen," Jour. Avier. Chem, Soc, Vol. 73, p. 4112, Sep- 
tember, 1951. 
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Earlier Work on '*Real" Surfaces 

While the above results are applicable to surfaces in contact with 
oxygen alone, it is obvious that etches commonly used to clean the 
germanium surface will lead to chemisorption of etch constituents so 
that we are dealing with a contaminated surface. On such a surface the 
contamination may significantly change the rate of further oxidation 
in a manner well known from oxidation studies*^ once the germanium 
has been removed from the etch. Furthermore, as most germanium 
etches are oxidizing, they will leave an initial thickness of oxide on 
the surface. It has been found by Green and Smythe^ that an etcherl 
surface may be covered by three monolayers of oxide five minutes after 
etching, increasing to six monolayers after three days' exposure to 
room air. Contaminated surfaces are presently less understood than 
clean surfaces, and the considerable amount of data that exists on 
surface traps has not yet been interpreted and correlated with th-e 
chemical and physical processes that are responsible for the traps. 
The present article constitutes a step in this direction. 

Two Basic Surface Treatments 

The results of two well-investigated surface treatments for gei- 
manium are taken as a starting point. These results are used to form 
a hypothesis about the oxidation of germanium which is then checked 
against the behavior of germanium transistors on life test. The two 
surface treatments referred to are the so-called No. 5 etch* and anodic 
oxidation in glacial acetic acid. The No. 5 etch is known to leave 
germanium monoxide, GeO, on the germanium surface, thereby making 
the surface conductivity strongly n-type.^ This gives p-n-p transistoi-s 
a low surface recombination velocity and consequently a high a„f, and 
a low saturation current. 

The anodic oxidation treatment in acetic acid has been worked oui 
by Zwerdling and Green'' and is known to leave a layer of germanium 
dioxide, GeOo, on the surface. This dioxide layer influences the surface 
conductivity from n-type towards less n-type as is shown in the foUow- 

^ R. K, Hart, *'The Oxidation of Aluminium in Dry and Humid Oxygen 
Atmo.sphere/' Froc, Roy Soc. A, Vol. 23G, p. G8, July,' 1956. 

^ M. Green and K. Smythe, as reported in Reference (1). 

* The No. 5 etch, developed by S. G. Ellis of these Laboratoi-ies, consists? 
of 40 parts 48 7r HF, G parts 30% H,0., 24 parts H^O. 

^ S. G. Ellis, private communi'^ation. 

9 S. Zwerdling and M. Green, ''Anodic Oxide Films on Germanium/' 
Joifv, Electrocheni. Soc, Vol. 103, p. 61, March, 1956 (abstract #60). Also 
Lincoln Lab. M.LT., Quarterly Progress Report, August 1, 1956. 
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ing experiment. Four groups of p-n-p transistors were etched electro- 
lytically in KOH and subsequently two of the groups were given the 
anodic oxidation treatment. Thereafter all groups were put on 85 °C 
shelf life test. Figure 1 is a plot of a^^^ versus time. Each experimental 
point is an average of three units. As can be seen in Figure 1 the 
untreated control groups behave quite normally, the experimental points 

100 1 1 




TRANStSTORS 



20 - 



0 10 20 30 

HOURS 

Fig, 1 — Influence on transistors of anodic oxidation in glacial acetic acid 
plus sodium acetate. In the ordinate scale, 100 represents the value of olci, 
immediately after the electrolytic etch. The solid curve represents Equa- 
tion (10). 

falling around the theoretical curve. The oxidized groups first show 
very low a^.^^ which, however, climbs rapidly (due to removal of bound 
etch constituents) and approaches a curve similar to the solid curve 
but shifted to the left. In other words, the oxidation treatment has 
applied an oxide layer of the same thickness as would be applied by 
heating to 85'' for some 30 hours. 
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The terms GeO and GeO^ have been used loosely as it is known 
that these oxides under the prevalent conditions take up a rather 
indefinite amount of H^O forming various hydrated forms such as 
di- or metagermanic acids, etc.^'^ For silicon, which behaves in a 
similar manner, these processes have been studied, but for germanium 
the details are not as well known. Therefore the terms monoxide (GeO) 
and dioxide (GeOj) have been used for simplicity without specifying 
amount of hydration or atomic arrangement. 

Without offering further proof at this time, it will be assumed 
that similar oxidation processes take place during the conventional 
processing of transistors. In particular, it will be assumed that the 
high a^i of p-n-p transistors found after an electrolytic etch in KOH 
is a result of a layer of GeO on the surface of the germanium, while 
the subsequent decay of cc,.,, may be interpreted as due to a growing 
layer of GeO^i on top of the GeO layer. In line with this it has been 
found by E. 0. Johnson and G. C. Dousmanis of these Laboratories 
that a germanium surface freshly etched electrolytically in KOH is 
n-type, but in the course of many hours gradually goes over more and 
more towards less n-type. 

That oxidation of the surface should influence the electrical char- 
acteristics is entirely in line with the present view of oxidation 
processes,^^'^- In the adsorption of gas atoms on the surface, strong 
electric fields are set up, as witnessed by the well-.known changes in 
surface recombination velocity, surface conductivity, etc., caused by 
exposure to different gas species.^ The strong electric field is funda- 
mental in the early growth of an oxide layer by aiding the transport 
of ions through the layer. As the layer grows, the electric field at the 
metal-oxide interface necessarily decreases, and it is this gradual 
reduction of the electric field that is responsible for the change in 
surface recombination velocity and thereby also a,.,, of transistors. 

Experimental Results 

The experimental work consists of a number of life tests on tran- 
sistors, comprising measurements of versus time under various 
conditions. 

0. H. Johnson, "Germanium and its Inorganic Compounds," Chem. 
Rev., Vol. 51, p. 431, December, 1952. 

11 K. Hauffe, ReaJdioveji in mid an Festen Stoffen, Springer Verlag, 
Berlin, Germany, 1955, p. 411. 

i-W. E. Garner, Chemistry of the Solid State, Butterworths Scientific 
Publications, London, England, 1955, p. 336. 
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0 40 80 120 160 200 240 280 

HOURS ^ 

Fig. 2 — Life test results for encapsulated silicone-covered p-n-p g-ermanium 
transistors. Shelf life at 70"C. (F. L. Hunter) 

It may be pointed out here that this data should not be interpreted 
as the performance of commercially available transistors. Factory 
procedures generally include an ''aging" step, e.g., heating to an ele- 
vated temperature for a period of time, which reduces a^, to ai)proxi- 
mately 35 per cent of its original value and after which moderate 
temperatures within the rating of the transistor do not cause signifi- 



160 9 




80 120 160 
HOURS ^ 

Fig, 3— Shelf life at 85''C. 
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HOURS ^ 

Fig. 4— Shelf life at 105"C. 

cant further changes. Thus commercial processes employ in sonf>e 
degree an oxidation process to achieve stability. 

Figures 2-5 present life test data for four groups of p-n-p ger- 
manium transistors at four different temperatures. In each group there 
is one transistor for which a^^^ seems to decay somewhat differently 
from the other three, perhaps because of a slight error in the initial 

reading. Therefore, an average has been taken of the three similar 
units, disregarding the fourth. In order to let all transistors contribute 
equally to the average, relative values have been used, setting the initial 

— The I'esults are shown in Figure 8. For each point, the 

160 : , , ^ , : , , : ^ . 1 ' 1 
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0 40 80 120 160 200 240 280 
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Fig. 5— Shelf life at 125"C. 
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spread between the three values is shown. The transistors had been 
electi^olytically etched in KOH and hermetically encapsulated with a 
silicone resin covering the germanium surface. 

Other runs were made in which it was established that the silicone 
was inconsequential and that units in silicone or air showed identical 
behavior within the accuracy of the measurements. Vacuum-sealed 
units showed a reduced decay as expected if the decay is caused by 
an oxidation process. Figure 6 represents such a run in which the 
units were sealed at lO--* mm Hg and life tested at lOO'^C, 70°C, and 
room temperature. Each point is an average of three units. Relative 
ralues are used as explained above. The theoretical curves are identical 
to Equation (10) except for a factor 0.09 representing the influence 
of reduced oxygen pressure. 

In order to test the theory describing these results, especially the 
temperature dependence of the oxide growth, a room-temperature life 
test was run ; this is shown in Figure 7. These transistors were sealed 
in hermetic cans in dry air without silicone. The results are also shown 
in Figure 9, where each point represents an average of five transistors, 
and again the spread is shown. In analyzing this test, the three tran- 
sistors with low have been discarded. 



The thickness of the growing GeOo layer cannot be very large at 
the temperatures (25° to 125°C) and times involved. It is shown 
later that this growth cannot exceed one hundred monolayers. It is 
reasonable, therefore, to assume that the increase of the surface recom- 
bintion velocity, s, is directly proportional to the increase of the 
thickness of the layer. For « u\ 



where iv is the thickness of the GeOo layer and v is a proportionality 
factor. Upon integration, 



is directly related to the surface recombination velocity by the 
equation^^ 



13 \V. M. Webster, "On the Variation of Junction-Transistor Current- 
Amplification Factor with Emitter Current," Proc, LR,E,, Vol. 42, p. 914, 
June, 1954. 



Outline of Method of Analysis 



ds dw 
dt ^ dt 



(4) 



S — So + V7V. 



(5) 
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Fig. 6 — Life test results for vacuum-sealed (10--* mm Kg) transistors. 
Shelf life at 105 C. (R. R. Johnson) 
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Fig. 7 — Life test results foi- transistors sealed in dry air, no silicone. 
Shelf life at room temperature. (E. L. Jordan) 
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Fig. 8 — Figures 2-5 replotted using average of relative values of a'ct.. 
Kange within which values fall is indicated. 




Fig. 9 — Figure 7 replotted using average of relative values of oLch. 
Range within which values fall is indicated. 
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1 

=C5S + Co, (6) 

where Cr, and C(; are given by the geometry, material, etc. of the tran- 
sistor. Combination of Equations (5) and (6) then gives 



= +CrJ)W, (7) 

where 

1 

'- "o 

Now by insertion of experimental values of a^.^, versus time in Equation 
(7), the growth rate of the oxide layer can be obtained from the 
resulting relation of w versus time. 



Analysis of the Data 

It is known that within a certain temperature range the oxidation 
process follows one of a small number of possible rate eiiuations.^^-^^ 
Using Equation (7) above, it is now possible to fit the rate equations 
to the experimental data. Tt, is then found that the so-called **cubic" 
rate equation, 

d = CjV^ + Cs. (S) 



where 7i — 1/8, fits the data very well while none of the other rate 
e(iuations can account for the initial drop. In Figures 10 and 11 
experimental points from Figure 8 have been redrawn and a number 
of possible I'ate equations have been fitted at two points. As can be 
seen, good fit can only be claimed for the **cubic." 

Equation (8) applies to the case of oxidation of a metal with a 
thin oxide layer. The validity of the assumption of a thin oxide layer 
(< 500 a) will be demonstrated by the dii'ect oxidation experiment 
described below. 

In Equation (8), the coefficient contains an exponential factor 
characterizing the activation energy for oxidation. Thus Equation (S) 
can be written as 



= C.:,e ^' t'A + Cr, (9) 

«c6 
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Fig. 10 — Some known oxidation equations fitted to the 70'' C 
curve of Figui'e 8. 
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Fig. 11 — Some known oxidation equations fitted to the 105°C 
curve of Figure 8. 
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Equation (9) has been fitted to the experimental points in Figure 8 
using coincidence at two arbitrary points, namely 20 hours and 160 
hours. The result is the equation 

=:95e ^' ^J^3 + 0.0091, (10) 

which is represented by the curves of Figure 8. As can be seen in 
Figure 9, this equation fits the room-temperature curve also. 

We can press the fit a little further. At a first glance, it is sur- 
prising that the theoretical curves, Equation (10), do not go through 
the point (0,100), in other words, that the constant in Equation (10) 
is not 0.01. However, it has been found by R. R. Johnson of the RCA 
Semiconductor Division, that when transistors are taken oflf a high- 
temperature life test at, say, 85°C and brought to room temperatur*? 
for the purpose of measurements, af.,, shows a decline for about 48 
hours as shown in Figure 12. This was not known at the time of the 
life test described above and therefore no time interval was fixed 
between removal from life test and measurement. The procedure was 
to remove the units in the morning and to measure them at any time 
during that day. It may therefore be assumed that an average of four 
hours elapsed between removal and measurement corresponding to a 
decline in a^.^, of about 9 per cent as shown in Figure 12. (This fact 
at the same time may explain some of the spread of the measured <y^,j, 
values.) However, the initial a^.,, was measured before any heat treat- 
ment and therefore would correspond to the end value in Figure 12'. 
It may be expected then that all the points except the initial reading 
are about 9 per cent high as indeed seems to be borne out by Equation 
(10). When the measurements are corrected for this effect they obe\ 
Equation (11) : 

3000 

~ tV^ + 0.01. (11) 

The value of the activation energy in Equation (10), (7200 cal/mol) 
at first seems unusually small. However, recent work on the oxidation 
of copper, which according to the reasoning outlined above follows 
a similar oxidation process, also suggests a small activation energy j, 
about 9,000 cal/moL 

" F. W. Young, J. V. Cathcart, and A. T. Gwathmey, ''The Rates of 
Oxidation of Several Faces of a Single Crystal of Copper as Determined 
with Elliptically Polarized Light," Acta MetaUurgica, Vol. 4, p. 145, No. 
2, 1956. 



1 : 

= 95e 

«c6 
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Direct Oxidation Experiment 

In order to provide a direct proof of the oxidation hypothesis set 
forth above to explain slow decay of of p-n-p transistors, the 
following experiment was performed. Samples of n-type germanium 
were heated in an oven at 85 °C for various periods of time. Subse- 
quently the samples were washed in an NaOH solution to remove any 
oxide formed. The amount of germanium in the solution was then 



100 » 



80 - 



EACH POINT AN AVERAGE OF 
18 TRANSISTORS 




40 



20 - 



0 10 20 30 40 50 

HOURS-*- 

Fig, 12— Off-life decline of a,^ ("48-hour efTect'O- (K. R. Johnson) 

determined quantitatively by a spectrographic method. Figure IS shows 
the results. The experimental points are in reasonable agreement with 
the cube-root-of-time law suggested earlier. In order to retain the 
oxide in easily soluble form it was found necessary to keep the samples 
at 100 per cent relative humidity during the heat treatment. 

The results in Figure 13 make it possible to deduce a quantitative 
expression for the oxidation rate. Assuming that at the low relative 
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humidity present in encapsulated transistors the oxidation rate is not 
too different from the rate at 100 per cent relative humidity, the follow- 
ing expression gives an upper limit for the oxide layer thickness in 
A at 85°C: 

2V ^ 21 + 1.3)1/3, (12) 

It should be pointed out that this oxidation rate is valid for an 
etched, and consequently hydrated, surface and therefore may be higher- 
than for a clean germanium surface. 

Influence of a Sequence of Temperatures 
From the picture outlined above some predictions may be made con- 



OXIDE 
LAYER THICKNESS 



A 




0 I 1 1 ■ — 1 1 t 

0 I 8 27 64 125 HOURS 

(TIME) 

Fig. 13 — Thickness of oxide layer versus time. 

cerning the influence of oxidation on units that are subjected to a 
sequence of temperatures. As an example, consider transistors that are 
first aged at room temperature for 24 hours. Following this they are 
heated for 24 hours at 105°C. With the help of Equation (10) we 
may now follow the change in n',.,, due to oxidation. After the first, 
24 hours of room temperature aging, a^f, is down to 85 per cent of 
its initial value. After the next 24 hours at 105"C, a,^, is down further 
to 80 per cent of the initial value. From then on at a temperature of 
85°C, acj>,ci can be predicted from the following time dependent ex- 
pression : 

1 

— 0.00427 (t + 160) H + 0.033. (13) 
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Alternatively, assuming that the value of a^j, at the start of the SS^'C 
shelf life test is 100 per cent, and that all measurements (except the 
initial one) are made an average of 4 hours after I'emoval from the 
oven, we have 

1 1 

= [0.00427 {t + 160)^3 + 0,033]. (14) 

«cvcf 3.3 X 1.09 

This relation is plotted in Figure 14. In a similar manner the results 





Fig. 14 — Theoretical curve and experimental values of c^cb 
for shelf life at 85' C. 



of any combination of temperature and aging time can be pi^edicted. 
Figure 14 also shows some experimental points from life tests of units 
that have been subjected to the above-mentioned sequence of tempera- 
tures. 

The first constant in Equation (13) will be different for different 
types of transistors and depends on the ratio C^/Co in Equation (6). 

It should be pointed out once again that factors other than oxida- 
tion also influence the change of a,.^, during life. However, interpreta- 
tion of other effects is difficult unless the influence of oxidation is 
first singled out. 
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Conclusions 

In view of the evidence presented, it appears likely that in ger- 
manium p-n-p transistors which have been electrolytically etched in 
KOH and subsequently hermetically sealed, a^^i is influenced by oxida- 
tion as follows: 

1. The conventional electrolytic etch in 40 per cent KOH leaves 
a composite hydrated oxide layer, GeO — GeO^ on the surface. 

2. The high a^.^ resulting from this etch is mainly a consequence 
of the GeO layer. 

3. The GeOo layer will grow with time and the decay of «c6 is a 
consequence of this growth. 

4. The value of a^.^ after storage of a p-n-p germanium transistor 
at a temperature T°K can be predicted from 

3000 

«cbo 

l = C^e + 

where o',,^^ is the initial value after electrolytic etch, and Ci is 
a constant which depends on the geometry of the transistor. 
For a 2N77 transistor, C^ = 95. refers to the oxide thickness 
at the start of the storage period and is the time necessary to 
grow this thickness at the temperature T, may be obtained 
from 

a* 100 

^cbo 

l = C,e T t,y^^ 

where a^bi is the value at the start of the storage period. 
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PULSE-FIRING TIME AND 
RECOVERY TIME OF THE 2D21 THYRATRON^^^ 

By 

John A. Olmstead and Melvin Roth 

KCA Kleclron Tube Division, 
Harrison, N, J. 

Suntmary — TJie time-dejjendent or dynamic characteristics of the 2D21 
are of considerable importance to the application engineer. In many cases, 
these characteristics determine the ability of the associated circuit to func- 
tion properly. Three of the most important dynmnic characteristics of the 
2D2I are pidse-firing tinic, ionization time, and grid-recovery time. 

Pulse-firing time is the time required after application of the grid ptdse 
for the ionization in the grid region to reach a value such that the grid 
loses control of the breakdoivn '})rocess. In the 2D21 this time is esse7itially 
eqnal to the ionization time (time required for anode conduction). Pulse 
firing time is a function of anode voltage, and grid-cathode voltage, the 
significant parameter being the grid-cathode voltage. When this voltage 
is positive, the grid draws appreciable current and proper circuit analysis 
must be applied to deterniine pulse source requirements. 

Recovery time is the period required for the grid to regain control after 
anode current flow has been interrupted by removal of anode voltage. It 
is primarily a function of anode current prior to extinction, anode voltage, 
and the grid— cathode voltage. The actual grid^cathode voltage must be 
determined by jyroper circuit analysis, with particular attention given to 
the ion grid current flowing during deionization. 

With proper circuitry the 2D21 can be ionized in 1 microsecond, and 
grid control can be recovered ivithin 100 microseconds. The data presented 
shoidd enable the application engineer to design the necessary circuitry 
with a minimum of trial and error. Although the data is applicable only 
to the 2D2I, the general concepts apply to other thyratrons. 



T 



Introduction 

's^HE 2D21 thyratron has found wide application in the industrial, 
nailitary, and computer fields because of its small size and desir- 
able characteristics as a low-impedance switch. A great many 
of these applications involve the use of the 2D21 in pulse-type circuits. 
Because adequate data on the dynamic characteristics of the 2D21 has 
not been available, designers have not been able to predict the end 
performance of such circuits, and, consequently, have been compelled 
to design them largely by trial-and-error methods. 



* Portions of this paper are based on work done for a thesis submitted 
to Newark College of Engineering by J. A. Olmstead in June, 1957 as part 
of the requirements for the degree of Master of Science in Electrical Engi- 
neering, 
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This paper discusses and presents quantitative data on the two most 
important dynamic characteristics of the 2D21 -pulse-firing time and 
grid-i'ecovery time. 

Pulse-firing time is the time required for the ionization in the gnd 
region to reach a value which causes the grid to lose control of the 
breakdown process. 

Grid-recovery time is the time required after removal of the anode 
voltage for the grid to regain a degree of control sufficient to prevent 
conduction when the anode voltage is reapplied. Although this internal 
is sometimes called the deionization time, the term recovery time is 
more precise because the grid can regain control even when a substan- 
tial amount of ionization is present. 

Pulse Firing 

'Tulse firing" of the 2D21 takes place in the following manner: 
Electrons emitted by the cathode are accelerated by the applied anode 
potential and/or grid potentials, and acquire sufficient energy to ionize 
the included gas. However, the region of the tube in which the electrons 
acquire ionization energies and the number which reach this region 
are determined primarily by the grid-cathode potential. At low 
grid-cathode potentials, electrons can acquire ionization energies only 
in the anode region, and few electrons reach this region. When the 
grid is made highly positive with respect to the cathode (for example, 
by application of a grid pulse) large numbers of electrons can acquire 
ionization energies in both the grid and anode regions. i 

If the grid pulse voltage produces sufficient positive ions in the grid 
region, these ions will surround the grid and neutralize the effects of 
the applied grid potential. As a result, the grid loses its ability to 
control the anode current and breakdown eventually occurs. 

The interval required for the grid to lose control after application 
of the firing pulse is called the pulse-firing time. The somewhat longer 
interval required for the tube to become completely ionized is called 
the ionization time. The pulse-firing time and ionization time of the 
2D21 are essentially equal and will be treated as such. 

Measurement of Pulse-Firing Time 

A block diagram of the circuit used to measure the pulse-firing time 
of 2D21, and the associated anode-voltage and grid-voltage waveforms 
are shown in Figure 1. The anode voltage is a 60-cycle square 



1 J. D. Cobine, Gaseous Conductors, p. 456, 457, McGraw-Hill Book 
Coinpany, Inc., New York, 1941. 
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TUBE FIRES 



Fig. 1- 



-(a) Circuit used to measure pulse-firing time of the 2D21; 
(b) and (c) anode and grid-voltage waveforms. 



wave. The applied grid voltage is a positive-going pulse delayed so 
that it occurs during the positive portion of this square wave. The 
adjustable series resistance R^i in the grid circuit is provided to facili- 
tate grid-current measurements. Figure 1(b) shows the grid and anode 
wave forms under conditions such that breakdown does not occur. 
Figure 1(c) shows the wave forms under breakdown conditions. 

Figure 2 is a time expansion of the breakdown waveforms in Figure 
1(c) showing the important parameters, namely peak grid-cathode 
voltage, E^;„ critical pulse width, and ionization time, T,-. 

The shaded portion of the grid pulse is unnecessary for firing the 



IONIZATION 
ANODE TIMECtO 




Fig. 2 — Expanded presentation of waveforms shown in Figure Ic, showing 
relation between ionization (pulse-firing) time and critical pulse width. 
Shaded portion of grid-voltage pulse is superfluous. 
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tube. In fact, the tube will fire even if the grid-pulse duration is 
slightly less than the ionization time. In the 2D21, however, the differ- 
ence between the critical pulse duration and ionization time is so sm^all 
that it can be neglected. Figure 3 shows the maximum I'equired 
grid-cathode voltage versus anode voltage for various ionization 
(pulse-firing) times. 

The curves in Figure 3(a) are for a shield-grid voltage of 0 volts 
and those in Figure 3(b) for a shield-grid voltage of H-9 volts. These 
curves can be used to determine the maximum grid-cathode voltage 
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Fig. 3 — Maximum grid-cathode voltage versus anode voltage with ionization 
time, Era, and Rg^ as parameters. 

required for any desired ionization time. For example, with 200 volts 
on the anode, and 0 volts on the shield grid, the grid-cathode voltage 
required to ionize the tube tvithin 5 microseconds is +17 volts. 

Note that the required grid-cathode voltage varies inversely w^ith 
ionization time and that the rate of variation is nonlinear. For example, 
the diflference between the grid-cathode voltages required for ioniza- 
tion times of 2 and 3 microseconds is substantially greater than the 
difference between those required for ionization times of 25 and 50 
microseconds. 
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Absolute time jitter due to variations in the grid-pulse amplitude 
will, therefore, be less at shorter ionization times. 

The grid-bias voltage may have any value within tube ratings; the 
important parameter is the actual grid-cathode voltage. 

Figure 4 shows curves of maximum grid current versus grid- 
cathode voltage, wath shield voltage as a parameter. These curves show 
the grid current flowing prior to firing. Any series impedance between 
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Fig. 4 — Maximum grid-current versus grid-cathode voltage, 
with Fjvc2 as a parameter. 



the pulse source and the grid, therefore, becomes an impoi'tant con- 
sideration in determining the actual grid-cathode voltage, because the 
applied pulse will be attenuated by the resulting voltage drop. The 
actual grid-cathode voltage for a given grid-circuit impedance and 
pulse-source voltage can be determined by circuit analysis. 

Figure 5 shows the minimum grid-cathode voltages required for 
various pulse-firing times and anode voltages. This information is 
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particularly valuable when the 2D21 is to be used in coincidence-type 
circuits because the curves can be used to determine the maximnm 
pulse width and/or amplitude that can be applied without firing the 
tube. For example, with 200 volts on the anode and 0 volts on the 
shield-grid, a pulse of 5 microseconds duration may fire the tube if 
the resulting grid-cathode voltage exceeds 8 volts. 
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Fig. 5 — Minimum grid-cathode voltage required for various ionization times, 
with anode voltage, Ecc^^ and Rg2 as parameters. 



Note that the data shown was obtained with the use of a resistance 
of 1,000 ohms in series with the shield grid. The use of this resistance 
reduces ion bombardment of the shield and thereby improves tube life. 
The value of 1,000 ohms was found satisfactory for use in all circuits. 



Recovery op^ Grid Control in the 2D21 

Once conduction has been established in the 2D21, flow of anode 
current cannot be interrupted or controlled by grid action. The grid 
loses control because it is surrounded by a sheath of positive ions 
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which neutralize the effects of any applied grid voltage. Although 
this sheath is relatively opaque to electrons, it does not fill the entire 
grid aperture, and, therefore, does not obstruct the flow of electrons 
to the anode. Increasing the negative grid voltage would normally 
expand this sheath and thereby permit the grid to control the anode 
current. During conduction, however, an increase in negative grid 
potential shifts the region of ion production in such a manner as to 
increase the ion density in the grid region rather than to allow sheath 
expansion. 

Ion production can be stopped only when the anode voltage is 
reduced below the extinction value (approximately zero volts). The 
resulting elimination of ion production permits the ion sheath sur- 
rounding the grid to expand. When the sheath expands suflftciently 
to close the grid aperture the grid begins to reassert control. If the 
grid potential at this point is sufficiently negative, subsequent reappli- 
cation of anode voltage will not refii^e the tube. As deionization pro- 
gresses, the grid acquires greater degrees of control, until it reaches 
its prefiring condition. Recovery time, therefore, is a function of the 
grid-cathode voltage at the instant anode potential is reapplied. This 
actual grid-cathode voltage is not necessarily the bias voltage. 

To determine the actual grid-cathode voltage during the deioniza- 
tion period, it is necessary to know the magnitude of the ion current 
flowing in the grid circuit throughout this period. This ion current 
(called the '*after-glow grid current*') depends upon the ion density 
within the tube and consequently, on the anode current prior to extinc- 
tion, and is, therefore, substantially independent of the external grid- 
circuit parameters. The after-glow grid current decays approximately 
exponentially from its initial value. If the voltage drop across the 
bias impedance due to the after-glow current tends to exceed the bias 
voltage, electrons as well as positive ions flow to the grid, cancelling 
a portion of the after-glow current and thereby preventing the grid 
from becoming positive with respect to the cathode. When the after- 
glow current decays to a value such that the resulting voltage drop 
across the bias impedance is equal to the bias voltage, the grid potential 
begins to return to the bias value. The resulting delay in the return 
of the actual grid-cathode voltage to the bias value is responsible for 
the increase in recovery time observed when a grid-bias impedance 
is present.^ 

- L, Malt^r and E. 0. Johnson, "Studies of Thyratron Behavior," RCA 
Review, Vol. XI, pp. 165-189, June, 1950. 
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Measurement of Recovery Time 

A block diagram of the setup used to measure the recovery time of 
the 2D21, and the associated waveforms are shovt^n in Figure 6. The 
test method shown is similar to one now under consideration by the 
Joint A.I.E.E.-I.R.E. Subcommittee on Gas Tubes. A master oscillator 
places repetitive positive pulses simultaneously on the control grid and 
anode. The tube conducts for a period sufficiently long to establish 
steady-state conditions. Recovery time is then measured by means of 
a "probe" pulse which may be applied to the anode at any desired time 
after cessation of the anode-current pulse. The delay of the probe 
pulse with respect to the end of the anode-current pulse is set at the 
desired recovery time. With the anode-conduction current and the 
probe-pulse voltage at the desired values, the grid-cathode voltage is 
then adjusted until the tube just fails to conduct during the probe pulse. 

Measurement of the after-glow current is facilitated by resistor 
Rg^ in the grid circuit. 

Figure 7 shows the critical grid-cathode voltages required for 
various recovery times. Figure 7a is for a peak anode current of 12 
amperes and 7b for a peak anode current of 5 amperes. Note that in 
Figure 7a after 120 microseconds of - deionization, a grid-cathode 
voltage of— 68 volts is required to prevent refiring at an anode potential 
of 200 volts. In Figure 7b, for the same anode voltage and recovery 
time, the grid-cathode voltage required to prevent refiring is now only 
—38 volts. 

Figure 8 shows similar curves for peak anode currents of 1 ampere 
and 0.5 ampere. Note that as peak anode current decreases, the required 
grid-cathode voltage for a given recovery time also decreases. It has 
been found that anode currents of less than 0.5 ampere have little effect 
on critical grid-cathode voltage. The curve for 0.5 ampere may, there- 
fore, be used for lower currents. 

Figure 9 shows the decay of after-glow current with time. These 
curves make it possible to determine the instantaneous ion grid current 
at any time during deionization. The actual grid-cathode voltage may 
be determined by circuit analysis. This analysis applies even though 
the resultant grid-cathode voltage sometimes appears to go positive. 

The magnitude of the positive-ion current flowing to the grid during 
conduction is the magnitude at zero deionization time. This current 
can be used to determine the grid-cathode voltage during the conduc- 
tion period; it should not exceed the — 10-volt rating for the 2D21. 
An example illustrating the use of the data presented is given below. 

The circuit under consideration is shown in Figure 10. The known 
parameters are: anode-supply voltage {E^^} =200 volts, grid-supply 
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Fig. 6 — (a) Circuit used to measure gTid-recovery time of the 2D21; 
(b) and (c) anode and grid-pulse Avaveforms. 
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Fig. 7 — Critical grid-cathode voltage as a function of recovery time, anode 
voltage, and Rg2, for peak anode currents of 12 and 5 amperes. 
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voltage {E^^{)—' 75 volts, peak anode current (1^=12 amperes, 
ionization time (T,,-) =5 microseconds max., grid-recovery time (T^) 
= 120 microseconds max. 

The problem is to determine the maximum value of Rg^, and the 
minimum pulse-source voltage required to fire the tube within 5 
microseconds. The value of Rg^ is determined from the data given in 
Figures 7a and 9. In Figure 7a, the intersection of the curve for 
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(b) 

Fig. 8 — Critical grid-cathode voltage as a function of recovery time, anode 
voltage, and Rg2j for peak anode currents of 1 and 0.5 ampere. 



T,. = 120 microseconds with the coordinate Ei,^ — 200 volts shows that 
a grid-cathode voltage of —68 volts at the instant E^,^ is reapplied will 
prevent firing. In Figure 9, the intersection of the curve for = 12 
amperes with the 120-microsecond coordinate shows that the after-glow 
current flowing in the grid circuit 120 microseconds after conduction 
ceases is 0.15 milliampere. 

From these values (^^^i^^-GS volts, and i,i = 0,15 milliampere) 
the maximum value for Rg^ is determined as follows: 
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Fig. 9 — ^'After-glow" grid-current as a function of recovery time 
and peak anode current. 
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Fig-. 10 — Basic circuit showing parameters involved in sample 
calculation of Rgi and E,^. 
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E. 



rrl 



75 — 68 



46,500 ohms. 



0.15 X 10-^ 



The next unknown to be determined is the minimum pulse-source 
voltage E^^ required to fire the tube within 5 microseconds. This voltage 
is found from the data shown in Figui'es 3a and 4. In Figure 3a, the 
intersection of the curve for = 5 microseconds with the E^^ = 200 
volts coordinate show^s that a grid-cathode voltage of +17 volts is 
required to fire the tube. In Figure 4, the intersection of the curve 
for Ey2 ~ 0 volts with the -f 17-volt coordinate show^s that the prefiring 
grid current is 0.29 milliampere. 

From these values {E^^^^ = +17 volts, i^^ = 0.29 milliampere, 
= 46,500 ohms, and E^^j — —75 volts) the pulse-source voltage 
is found as follows: 



When the 2D21 is to be used in applications requiring short ioniza^ 
tion and recovery times, the external ciixuit parameters must be care- 
fully chosen. The important parameter is the grid-cathode voltage. 
To determine this voltage, knowledge of the grid currents, and careful 
analysis of the external grid circuits are necessary. Generally grid 
impedances should be kept at a minimum. Large grid signals and bias 
voltage may be required. 

Ionization times as short as 1 microsecond may be achieved in the 
2D21 by pulsing the grid sufficiently positive. Because appreciable grid 
current flows under these conditions, the effective impedance of the 
signal source becomes a very important consideration. This impedance 
attenuates the pulse and should, therefore, be kept as low as possible, 
preferably below 50,000 ohms. 

Grid control in the 2D21 can be recovered in about 100 microseconds 
if the circuit is properly designed. The important parameters affecting 
recovery time are anode current prior to extinction, anode voltage, and 
the grid-cathode voltage. The grid-cathode voltage is determined by 



E,,j,^i,,R^,+E, 



= 17 + (0.29 X 10-3 X 46,500) + 75 106 volts. 



Conclusion 
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the after-glow grid current, grid-circuit impedance, and bias voltage. 
A grid-circuit impedance below 50,000 ohms and large bias supply 
voltages are recommended. 

The data presented covers the ranges of operating conditions nor- 
mally employed for the 2D21, and should enable the design engineer 
to predict the end performance of his circuit with reasonable accuracy. 
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Correction: 

In the paper entitled "Validity of Traveling-Wave-Tube Noise 
Theoiy," by W. R. Beam and R. C. Knechtli, which appeared on pages 
24-38 of the March 1957 issue, the following information was inad- 
vertently omitted : 

The authors wish to thank R. W. Peter for his guidance. 
Part of the work described in this paper was performed under 
a Signal Corps Contract. 
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Jack Sklansky received the B.E.E. degree from the City 
College of New York in 1950, the M.S.E.E. degree from 
Pur due University in 1952, and the Doctor of Engineer- 
ing Science degree from Columbia University in 1955. 
While at Purdue he served as a teaching assistant. At 
Columbia he was a Eugene Higgins fellow from 1952 to 
1954 and a research assistant at the Electronics Research 
Laboratory during two summers and 1954-55. During 
previous summers he was associated with the U. S. Bu- 
reau of Reclamation, the U. S. Naval Ordnance Labora- 
tory, and the Bell Telephone Laboratories. Since 1955, 
he has been a member of the technical staff of the David Sarnolf Research 
Center, where he has been investigating problems in satellite instrumenta- 
tion, missile dynamics, and radar-data processing. Dr. Sklansky is a member 
of the Institute of Radio Engineers, the American Institute of Electrical 
Engineers, Tau Beta Pi, Sigma Xi, and Eta Kappa Nu. 



J. ToRKEL Wallmark received the Civilingenjor degree 
in Electiica! Engineering from Royal Institute of Tech- 
nology, Stockholm, Sweden, in 1944, the Teknologie 
Licentiat degree in 1947 and the Teknologie Doktor 
degree in 1953. He was a research assistant with the 
Institute from 1945-1953 and also held positions with 
various other organizations active in research in electron 
tube and semiconductor problems. From 1947 to 1948 
Dr. Wallmark was a fellow of the American- Scandi- 
navian Foundation at RCA Laboratories, to which he 
returned in 1 953 ; he is presently w^orking on semi- 
conductor surface problems. 
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